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ABSTRACT 


This  report  describes  in  detail  the  use  oi  the  Project  G1  7  radar 
lor  making  incoherent  scatter  measurements  ol  the  ionosphere.  The  data 
processing  and  analysis  methods  are  lully  explained,  including  all  the 
assumptions,  interpolations,  and  extrapolations.  'IVo  new  major  computer 
programs  were  implemented.  Portions  ol  the  radar  hardware  that  are  es¬ 
pecially  important  to  the  incoherent  scatter  measurements  and  portions 
that  have  been  developed,  modilied,  or  improved  during  the1  past  year  are 
also  described.  The  development  and  implementation  ol  a  digital  auto¬ 
correlator,  in  combination  with  new  computer  soltwarc,  has  resulted  in 
an  order-ol -magnitude  improvement  in  data  processing  time,  as  well  as  a 
significant  improvement  in  measurement  accuracy. 
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I  INTRODUCTION 


The  incoherent  scatter  technique  has  made  possible  the  measurement 
of  ionospheric  electron  densities  and  electron  and  ion  temperatures  with 
a  single,  g i ound-based  radar  system.  These  ionospheric  parameters  are 
important  to  the  Defense  Atomic  Support  Agency  (UASA)  nuclear  weapons 
test  program.  During  the  past  live  years,  Stanford  Research  Institute 
has--for  DASA--designed ,  constructed  and  operated  an  incoherent  scatter 
radar  to  develop  the  techniques  for  operating  in  a  nuclear  environment, 
to  improve  the  radar  system,  and  to  conduct  ionospheric  researches  on 
the  natural  ionosphere. 

During  the  past  twelve  months,  the  objectives  of  this  project  have 
been  to  implement  a  program  of  equipmental  and  operational  improvements 
to  the  Mighty  Sky  Project  617  radar  in  order  to  achieve  a  high  degree  of 
test  readiness.  Specifically,  the  project  efforts  during  the  last  half 
of  1969  and  the  first  quarter  of  1970  have  been  directed  towards: 

(1)  Implementing  computer  programs  and  procec ures  to  speed 
up  and  automate  data  processing  and  data  analysis 

( 2 )  Implementing  a  digital  autoco  relator  to  enable  spectral 
analysis  to  be  performed  in  (nearly)  real  time 

(3)  Conducting  ionospheric  studies  using  the  incoherent 
scatter  technique  for  the  purpose  of  better  under¬ 
standing  ionospheric  phenomena 

( 4 )  Planning  the  move  of  the  radar  to  the  auroral  zone 

(b)  Maintaining  a  current  awareness  in  the  progress  of 

the  incoherent  scatter  technique  so  that  maximum  use¬ 
fulness  of  this  technique  for  DASA's  purposes  will  be 
realized , 


'Jliis  report,  Final  Report --Pa i*t  A,  describes  the  l-csults  of  the 
work  perfonned  to  satisfy  Objectives  1  and  2.  It  is  also  intended  for 
use  as  a  reference  document  t ha t  describes  in  detail  the  hardware, 
software,  and  procedures  used  in  processing  and  analyzing  the  incoherent 
scatter  data.  Since  many  changes  in  hardware  and  software  have  been 
made  during  the  past  five  vears.  and  since  the  radar  will  be  relocated 
to  Alaska  and  used  by  personnel  less  familiar  with  its  operation,  it  was 
felt  that  a  report  collecting  in  one  volume  all  the  pertinent  information 
concerning  the  receiver  system,  calibration  system,  analog- 1 o-d igi t a  1 
conversion  (ADc)  system,  digital  subsystem,  and  computer  processing 
system  would  be  extremely  useful. 

The  work  performed  under  Objective  3  of  this  project  is  described 
in  a  separate  report,  Final  Report — Part  B.  That  report  includes  the 
results  of  three  48-hour  data  runs  (June,  July,  and  August  19G9)  taken 
before  the  installation  of  the  digital  autocorrelator  and  at  least  three 
48-hour  data  runs  (February,  March,  and  April  1970)  taken  after  the  in¬ 
stallation  of  the  autocorrelator.  In  addition,  data  were  taken  and  anal¬ 
yzed  for  two  partial  solar  eclipses — one  on  11  September  19G9  in  which 
the  sun  was  about  70  percent  obscured,  and  one  on  7  March  1970,  in  which 
the  sun  was  less  than  30  percent  obscured. 

The  work  perfonned  under  Objective  4  of  this  project  has  been 
completed;  it  resulted  in  a  proposal  entitled  "DASA  Project  G17  Radar 
Relocation,"  which  was  submitted  to  DASA  in  early  March  1970.  Included 
in  that  proposal  is  a  description  of  the  work  done  under  Objective  4. 

Objective  5--maintaining  a  current  awareness  of  the  progress  of  the 
incoherent  scatter  technique — is  applicable  to  all  the  other  objectives 
of  this  project.  The  achievement  of  this  objective  is  evidenced  by  the 
improvements  and  results  described  in  botli  Part  A  and  Part  B  of  this 
Final  Report. 
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The  plasma  physics  theory  governing  incoherent  scattering  of  radio 

1  2  3  * 

waves  from  the  ionosphere  is  well  described  in  the  literature  *  ’ 
and  will  not  be  repeated  here.  However,  the  equations,  assumptions, 
interpolations,  and  extrapolations  used  in  the  data  processing  and 
analysis  are  described  in  Section  II,  Measurement  Technique.  That 
section  gives  the  mathematical  procedures  used  in  the  processing  of 
the  incoherent  scatter  data  so  as  to  yield  electron  densities  and 
electron  and  ion  temperatures. 

Section  111  of  this  report  describes  the  radar  system  hardware, 
inc luding : 

(1)  Receiving  and  calibration  equipment 

(2)  Analog- to-d igi tal  conversion  equipment 

(3)  Digital  autocorrelator 

(4)  Computer  interfacing  equipment. 

The  software  developed  for  and  used  by  the  Project  (517  radar  is 
described  in  Section  IV.  Included  are  the  on-line  data  collection 
program  and  the  off-line  data  processing  and  analysis  program,  both 
of  which  were  implemented  during  the  past  nine  months. 

Section  V  summarizes  the  improvements  incorporated  into  the  system 
during  the  past  year  and  points  out  further  improvements  that  could  be 
made  in  the  future. 


*  References  are  listed  at  the  end  of  this  report. 
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11  MEASUREMENT  TECHNIQUE 

Incoherent  scatter  is  basically  a  technique  l'or  measuring  ionospheric 
electron  densities  and  electron  and  ion  temperatures  by  using  a  ground- 
based  radar  facility.  Several  alternative  measurement  methods  are  pos- 

4 

sible.  This  report  discusses  in  detail  only  the  monostatic  pulsed  radar 
technique  implemented  for  the  DASA  Project  G17  L-band  radar.  Measurements 
of  received  power  as  a  function  of  range  and  ol  the  frequency  spectrum  of 
the  received  signal  as  a  function  of  range  are  used  to  determine  the  iono¬ 
spheric  parameters. 


A.  Power  Profile 

Equation  (l)  is  the  formulation  of  the  radar  equation  applicable  to 

c; 

the  incoherent  scatter  measurements:' 


where* 

/  d. 

N  =  electron  density  at  range  K  (el/m  ) 

c  =  velocity  of  light  (m/s) 

-28  2 

0  -  0.5  X  10  m  =  backscntter  cross  section  ol  single 

m  ■> 

electron,  assuming  Q“  -  0,  T  =  1 

r 

All  terms  are  listed  and  defined  in  the  Glossary,  which  precedes  the 
Kef erences . 
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G  =  antenna  gain  along  main  axis  =  47.15  dil 
o 

tG  —  antenna  halt  beaniwidth  (radians)  =  0.3  tt/180 
^  =  radar  wavelength  (meters) 

T]  =  product  oi  transmit  and  receive  system  losses  -  -1.9  dll 
1*  =  transmitted  power  (watts) 

T  -  transmitted  pulse  length  (seconds) 


P 

K 


received  power  at  range  1’  (watts) 


T  =  electron-to-ion  temperature  ratio  at  range  R 
r 

•>  v  (5 

a  =  (kl>)“  =  14.22  X  10  '1  /N  at  X  =  0.23  m. 

e 

'liie  nunierical  constants  (first  bracket)  and  unchanging  radar  parameters 
( second  bracket)  can  be  evaluated  and  combined  to  give: 


N  = 


2  2 
( 1  +  cr  +'J  )  ( 1  +  cr  ) 

(o-  (^)  (»\) - r1 - 


In  this  form,  the  equation  is  separated  into  a  numerical  constant 
(first  bracket),  variables  that  may  be  changed  but  are  fixed  for  any  par¬ 
ticular  data  run  (second  bracket),  variables  that  arc  directly  measured 
(third  bracket),  and  variables  that  must  bo  inferred  from  other--f rcquency 

spectrum — measurements  ( fourth  bracket).  An  additional  complication  is 

v 

introduced  by  the  fact  that  cv  is  itself  a  function  of  N.  Thus,  lo  deter¬ 
mine  the  electron  density,  one  must  first  know,  or  at  least  estimate,  the 
electron  density.  The  solution  of  this  problem  is  discussed  later.  As  a 
first  stop  in  the  data  analysis  procedure,  the  density  is  computed  by  as sum- 
ing  a  =0  and  T  =  1,  giving  what  we  call  the  raw  electron  density. 


G 


/ 


(3) 


N 


w’'erc 

20 

C  -  0.475  X  10  =  .system  constant, 

s 

To  determine  the  electron  density  from  Kq .  (2)  or  (3)  above,  a 

careful  measurement  of  the  magnitude  of  the  received  power  is  required. 
Also,  since  the  received  signal  is  at  best  only  a  few  times  larger  than 
the  noise  level,  and  at  most  ranges  loss  than  the  noise  level,  an  accu¬ 
rate  subtraction  of  the  receiver  noise  level  from  the  receiver  output 
containing  both  incoherent  scatter  signal  and  noise  must  be  made.  Further¬ 
more,  the  receiving  system  gain  changes  as  a  function  of  temperature  and  is 
periodically  adjusted  toy  the  operator.  The  received  noise  power  also 
changes  due  to  (l)  reception  of  signals  from  interfering  radars, 

(2)  rotation  of  the  celestial  sphere  through  the  antenna  beam,  and 

(d)  variations  in  the  gain  of  the  receiving  system  parametric  amplifiers. 


To  account  for  gain  and  noise  power  changes,  the  received  noise 
level  is  continuously  monitored  and  recorded.  The  noise  level  and  the 
signal-plus-noise  level  are  calibrated  by  periodically  injecting  a  noise 
pulse  of  accurately  known  amp  .itude  into  the  receiver  chain  in  front  of 
the  first  parametric  amplifier.  This  pulse  is  called  the  "calibrate 
pulse,"  Figure  1  shows  how  the  detected  received  output  as  a  function 
of  time  would  look  on  an  idealized  A-scope  display  if  a  square-law  de¬ 
tector  were  employed.  The  ordinate  is  average  received  power,  and  the 
abscissa  is  time  within  one  interpulse  period.  The  following  equation 

is  used  to  determine  the  received  power,  P  ,  at  range  R; 

r 


(4) 
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FIGURE  1  RECEIVED  POWER  AS  A  FUNCTION  01  TIME  (RANGE) 


whore 


I’  -  the  received  power  measured  at  the  j  range 
J 


P  =  the  received  power  measured  at  a  range  where  no  signal 
n 

is  present 


P  =  the  received  powei*  measured  at  ranges  where  the  calibrate 
c 

pulse  is  present 


P  =  the  injected  calibrate  noise  pulse  power  =  KT  1> 
k  k  1) 

B  =  the  receiver  noise  bandwidth  of  the  detected  channel. 

]) 

Bars  over  the  quantities  are  used  to  indicate  time  average.*..  Since  the 
incoherent  scatter  signal  is  noiselike,  all  measurements  must  be  done 
statistically,  i.e.,  quantities  must  be  measured  and  time  averaged  in 
order  to  determine  their  magnitudes.  Typically,  a  noise  pulse  tempera¬ 
ture  (T  )  oi  7 O'  K  is  used  i'or  our  measurements, 
k 
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The  received  powers  ( 1’  are  measured  digitally  by  sampling 

j  n  c 

with  an  analog-to-digi tal  converter  (ADC)  the  linearly  detected  received 
signal  voltage.  Squaring  of  the  voltages  to  obtain  power  and  averaging 
ol  the  power  is  done  digitally  in  real  time  in  a  general  purpose  digital 
compu ter. 

At  this  point,  a  comment  on  the  measurement  oi  radar  range,  H,  is 
in  order.  lor  a  hard  radar  target,  and  at  a  radar  lrequency  well  above 
the  plasma  lrequency  ol  the  ionosphere,  where  the  group  velocity  is  very 
close  to  the  velocity  oi  light,  the  range  can  be  measured  by  simplj  mean 
tiring  the  time  delay  between  the  leading  edge  ol  the  transmitted  pulse 
and  the  leading  edge  ol  the  target.  lor  this  situation, 

c  t  . 

K  =  ~  ,  (5 


where  t  is  time  delay,  and  the  factor  ol  2  accounts  lor  the  two-way 
propagation  path.  Tor  a  "soli"  (distributed)  target  such  as  the  iono¬ 
sphere  or  an  ionized  cloud,  the  returned  energy  represents  the  convolu¬ 
tion  oi  the  transmitted  pulse  with  the  target.  A  simple  calculation  ior 
this  case  shows  that  the  range  to  the  center  ol  a  distributed  target 
symmetric  in  range  is  given  by  the  time  delay  1 rom  the  center  oi  the 
transmitted  pulse  to  the  center  ol  the  target.  Thus  Tq.  (5)  must  be 
modified  for  incoherent  sc.it  ter  measurements  to  be 


where  T  i;-  the  transmitted  pulse  length. 


V.'e  can  now  summarize  the*  power  profile  measurements  as  follows. 

The  received  signal  is  linearly  detected  and  sampled  by  an  ADC  at  uni¬ 
form  intervals  (30  or  GO  Us  for  our  system)  during  the  interpulse  period. 
The  sampled  voltages  art;  fed  into  a  computer  where  they  are  squared  to 
give  the  P.'s;  j  is  the  index  indicating  the  radar  range  at  which  the 
sample  was  taken.  For  each  j,  the  l’.'s  from  many  consecutive  interpulse 

periods  are  averaged  to  give  P  's.  Then  P  and  P  are  calculated  by 

j  c  n 

averaging  the  P.'s  over  the  range  intervals  whore  the  calibrate  pulse 
is  injected  and  over  the  intervals  where  no  incoherent  scatter  signal 
exists,  respectively.  Next,  the  P  's  for  each  j  are  computed,  according 
to  Ft).  (*l).  From  the  P  's  and  K's  calculated  as  indicated  by  Kq .  (g), 
the  raw  density  N  is  calculated  by  using  Fq.  (3).  Then,  when  the  tem¬ 
peratures,  T  and  T  have  been  determined  (discussed  in  the  next  section), 
e  i 

the  true  density  can  be  calculated  by 

X  =  x'(l  -t  T  +  g2)(1  +  a~  )/2  .  (7) 

The  terms  in  parentheses  indicate  how  the  scattering  cross  section  changes 
v 

with  changes  in  aT  and  T  . 

r 

'this  met  hoi-  assumes  that  one  lias  very  accurate  knowledge  oi  the  radar 

system  parameters  that  go  into  the  system  constant,  C  ,  and  that  these 

parameters  do  not  change  or  arc  continuously  measurable  during  the  course 

of  a  data  run  and  from  day  to  day  and  month  to  month.  In  practice,  we 

have  found  it  is  easier  and  more  accurate  to  obtain  the  shape  of  the 

density  profile  as  discussed  above  but  to  normalize  the  maximum  density 

to  the  density  indicated  bv  measuring  f  F2  with  an  ionosonde. 

o 
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li. 


'1  empera  t  ure  Measurement  s 


1.  background 

'Jho  spectrum  oi  the  incoherent  backseat  ter  signal  is  ol  great 
importance  because  it  contains  all  the  inlormation  on  electron  and  ion 
temperatures.  'these  temperatures  are  of  interest  in  themselves,  and 
they  arc  also  necessary  for  proper  computation  oi  electron  density  pro¬ 
files  [hq.  (2)]. 

Figure  2  shows  the  spectrum  (that  is,  the  returned  power  per 
unit  bandwidth  plotted  as  a  function  of  Doppler  shift)  of  the  incoherent 
scatter  signal  for  a  transmitted  frequency  of  1290  MHz.  The  hall  spec¬ 
trum  is  shown  here,  with  the  total  spectrum  being  symmetric  around  zero 
Doppler  shift  in  the  absence  of  ionospheric  drills  or  currents.  Also 
shown  in  this  figure  arc  the  effects  on  the  spectrum  shape  of  various 
clcctron-to-ion  temperature  ratios,  T  .  As  T  increases,  the  "wing" 
of  the  spectrum  becomes  more  pronounced.  Therefore,  a  practical  way  of 
estimating  T  is  by  measuring  the  ratio  of  power  (wing)/powcr  (center 
frequency ) . 


0  4  8 

FREQUENCY - k Hz 


FIGURE  2  CALCOLATlD  INCOHERENT-SCATTER  IONIC  SPECTRA 
kD  0,  04  IONS,  T,  =  825  K,  Te/T,  VARYING 
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"111,.  j,in  1  L'lujH -r.i t arc  can  also  Ijc  deduced  ilom  the  measured  sp*c»,ra 

■  -1  - 

•flic  total  spectral  width  is  approximately  pruixirlion.il  to  [d^  ■*  1  }  *  -  - 
f  it;ure  3  slm.v;.  the  elleet  on  the  spectrum  ol  increasing  "I  while  keeping 
1  iixetl.  Thus  by  measuring  two  parameters  ol  the  spectrum,  say  peak-to- 

r 

valley  ratio  and  width,  one  can  theoretically  solve  lor  the  two  quantities 
ol  interes  t  ,  'I  and  1  .  . 


riCUfU  3  CALCULA1  LD  INC0III.RCN1  -SCAT  1  LR  IONIC 
M’LCIRA  kD  0.  0‘  IONS,  1,,/1,  =  2  0,  1, 
VARYING 


In  practice,  and  particularly  at  our  23-cm  wavelength,  an 

additional  complication  arises.  the  ratio  (o)  ol  the  bebye  length, 

[p  =  t) ( T  /N  )  '  to  the  operating  wavelength  allects  the  shape  oi 

e 

the  spectrum  as  well  as  the  backseat  ter  cross  section,  as  discussed 
earlier.  Figure  -1  shows  the  elleet  on  the  spectrum  oi  changing  a  b\ 
varying  X  but  leaving  the  ion  and  electron  temperatures  unchanged. 
Notice  that  both  the  bandwidth  and  the  peak— to— val ley  latio  ol  the 


spectrum 


are  a  Hoc  ted.  The  elleet  on  the  scattering  cross  section  can 


I 

<J 


— -  a  -  0.15  (N  -  I0  ‘  el/m  ) 


EIGURE  4  CALCULATED  INCOHERENT-SCATTER  IONIC 
SPECTRA  0‘  IONS  T,  =  800  K,  Tt,/T(  -  2.  kD 
VARYING 


At  this  point  wo  can  summarize  by  listing  the  three  physical 
quantities  that  aiiect  the  shape  oi  the  incoherent  scatter  spectra: 

(l)  electron  temperature;  (2)  ion  temperature;  and  (d)  electron  density 
11  one  oi  these  quantities  is  known,  the  other  two  can  be  deduced,  but 
because  oi  their  complex  interrelationships,  all  three  cannot  bo  ex¬ 
perimentally  deduced  lrom  spectrum  measurements  alone.  l'or  example, 
the  change  in  the  spectrum  shape  due  to  changing  N  in  Figure  4  can  be 


LiUhvo.)  uuju  i  on  l  t'u  u>  v-iiit  i 


Wii  lie  I  Oa  V  I  N 


i\ed  . 


2.  Instrumental  Liiccts 

Proper  measurement  and  analysis  oi  the  spectra  is  necessary 
lor  correct  interpretation  ol  ionospheric  densities  end  temperatures. 
Two  principal  iactors  complicate  the  measurement:  (l,  the  signal-to- 
noise  ratio  is  low--very  ol'ten  less  than  unity;  anil  (2)  the  measuring 


system  itself  (transmitter,  receiver,  and  digital  processor)  distorts 
tlie  spectrum.  Low  signal-to-noiso  ratios  limit  the  accuracy  and  height 
regime  over  which  measurements  can  be  made  and  lead  to  the  use  of  rela¬ 
tively  long  integration  times  (minutes  to  tens  of  minutes).  The  effects 
of  the  measuring  system  on  the  spectrum  can  bo  calculated  and  then  cor¬ 
rected  for  in  processing  the  data.  However,  the  measurement  accuracy 
is  somewhat  reduced,  since,  as  will  be  pointed  out  later,  the  equip¬ 
ments!  effects  are  such  as  to  make  the  inferred  parameters  T  and  T 

e  i 

more  sensitive  to  errors  in  spectral  shape  determination. 

Let  us  now  briefly  discuss  the  effects  of  the  measuring  instru¬ 
ment  (a  pulsed  radar  system)  on  the  measurement.  In  the  discussion  we 
make  use  of  the  fact  that  multiplication  of  two  functions  in  the  time 
(frequency)  domain  is  equivalent  to  convolution  of  their  Fourier  trans¬ 
form:.  in  the  frequency  (time)  domain.  Let  us  represent  the  power  spec¬ 
trum  of  density  fluctuations  in  the  ionospheric  plasma  by  .s(l'),  where 
f  is  the  Doppler  shift  frequency.  This  is  the  incoherent  scatter  spec¬ 
trum  that  would  result  If  the  measuring  equipment  did  not  perturb  the 
signal.  L'J  lie  equal  ions  for  calculating  s(  J  )  have  been  given  in  a  pre¬ 
vious  report,  and  Figures  2,  d,  and  d  are  some  representative  plots 
of  S(f )  vs.  I.]  'i  lie  quantity  s(  i  )  is  the  spectrum  that  would  result  from 
the  scattering  by  the  plasma  oJ  a  spectrally  pure  CW  signal. 

Since  :>  puls  >d  radar  is  used  lor  our  measurements,  which  means 
v/e  transmit  not  a  single  frequency  but  a  band  of  frequencies ,  the  fluc¬ 
tuation  spectrum  S( J)  inns t  be  convolved  with  the  power  spectrum  oi  the 
transmitted  pulse,  'the  result  oJ  this  convolution  operation  represents 
the  spectrum  oJ  the  incoherent  scatter  signal  at  the  receiving  antenna 
terminals,  kepresen t i ug  the  power  spectrum  oi  the  transmitted  pulse  by 
1  ’(  J  ) ,  the  signal  ,  S  f  J  ) ,  present,  at  the  antenna  terminals  is 

-\(J)  =  S(  J  )  *  )’(  J  )  ,  (8) 

where  the  asterisk  indicates  the  convolution  operation. 
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The  signal  S_^(  i )  then  passes  through  the  receiver  chain,  includ¬ 
ing  several  amplifiers,  mixers,  and  filters  which  have  bandpass  character¬ 
istics  that  affect  the  signal,  lotting  the  frequency  response  curve  of 
the  receiver  chain  bo  represented  by  I?(f),  the  signal,  S  (f),  after 
passing  through  the  receivers,  is  related  to  S  (f)  by 


S2(f)  =  8^1)11(1) 


(9: 


Next  the  signal  is  range-gated  to  obtain  spectra  at  a  number  of 
radar  ranges  (ionospheric  heights).  The  range-gating  process  can  he 
represented  in  tiie  time  domain  by  multiplication  of  the  signal  by  the 
gating  waveform  (a  square  pulse).  This  is  equivalent  in  the  frequency 
domain  to  convolving  S  ( f )  with  the  Fourier  transform,  G(f),  oi  the 
autocorrelation  function  of  the  square  pulse.  This  operation  gives 


S3(f)  -  S2(f)  *  G(f) 


( 10 ; 


where  S^(f)  is  the  incoherent  scatter  spectrum,  perturbed  by  the  measur¬ 
ing  instrument,  which  is  measured  by  the  radar  system. 

The  sequence  of  operations  described  above  is  indicated  graphi¬ 
cally  m  Figure  5.  Equivalent  operations  in  both  the  time  (autocorrela¬ 
tion)  domain  and  the  frequency  (power  spectrum)  domain  are  shown.  Pass¬ 
ing  from  one  side  to  the  other  of  the  dashed  vertical  line  indicates 
using  the  Fourier  transform. 

Figure  6  is  an  example  of  the  effect  on  the  spectrum  of  a  finite 

pulse  length  and  a  receiver  range  gate.  For  this  figure  the  receiver 

bandpass  has  been  assumed  to  be  flat.  Notice  that  the  greatest  change 

is  in  the  peak-to-valloy  ratio,  which  is  closely  related  to  T  /T  .  Also 

e  i 


FOURIER 


flGURE  5  SPECTRUM-PERTURBATION  BLOCK  DIAGRAM 
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FREQUENCY  -  kHz 

FIGURE  G  CALCULATED  INCOHERENT-SCATTER  SPECTRA,  SHOWING  EFFECT 
OF  FINITE  PULSE  W IDTH  AND  RECEIVER  RANGE  GATING 


changes  in  the  peak-to-val lev  ratio  for  a  given  change  in  T  /T  arc  not 

e  i 

as  great,  leading  to  reduced  ability  to  detect  small  changes  in  !  /'J' 

e  i 

and  thus  somewhat  reduced  accuracy.  It  the  spectral  perturbations  dis¬ 
cussed  in  this  section  were  not  allowed  ior  in  the  data  processing  and 

analysis,  one  would  underestimate  T  /'l'  and  overestimate  T  . 

e  i  i 


3.  Interpretation  tor  ct  =  0 


As  mentioned  earlier,  the  parameters  T  ,  T  ,  and  N,  as  well  as 

e  i 

the  pulse  length,  T,  and  the  range  gate  width,  w,  all  affect  the  shape 

oi  the  measured  spectra,  S  (f).  Throughout  a  data  run  T  and  w  remain 

fixed  and  can  be  accounted  for  in  the'  computation  of  the  theoretical 

spectral  shapes  used  for  comparison  with  measurement.  bet  us  ior  the 

2 

moment  assume  that  N  is  very  large  such  that  0t  ~  0  and  see  how  one  then 
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Later  we  will  show  haw  smaller 


interprets  the  data,  to  obtain  T  and  T  . 

c  i 

v 

densities  (larger  a  )  are  accounted  for  in  the  interpretation  of  the 
measurements . 

One  way  o::  obtaining  T  and  T  from  a  measured  spectrum  would 

e  i 

be  to  precompute  and  store  in  a  library  in  a  computer  a  vast  number  of 

theoretical  spectra,  for  many  values  of  T  and  1  .  The  spectra  so  com- 

o  i 

puted  and  stored  would  include  the  perturbations  due  to  the  pulse  length 

and  range  gate.  Then  one  could  have  the  computer  compare  the  measured 

spectrum  with  the  library  of  computed  spectra  and,  through  some  best-fit 

procedure,  decide  which  combination  of  T  and  T  best  matched  the  data. 

e  i 

This  technique  would  require  a  considerably  larger  computer  than  wo  have 
available  with  our  radar  system. 

An  alternative  method,  the  one  actually  used,  requires  much 
less  computer  time  and  storage.  it  involves  scaling  two  easily  identi¬ 
fiable  features,  the  peak-to-vallcy  ratio  (PV)  and  the  bandwidth  (liW), 
from  the  measured  spectra,  and  using  what  wo  call  a  spectrum  scaling- 

chart  to  interpret  these  parameters  in  terms  of  T  and  T  .  This  can  be 

r  i 

done  easily  and  quickly  in  our  computer.  Figures  7  through  10  show  the 
spectrum  scaling  charts  used  for  our  system.  They  were  prepared  for  the 
four  most  commonly  used  combinations  of  pulse  width  and  range  gate  width: 

T  =  360  ps  and  w  -  2-10,  TOO,  020,  400  Ps ,  From  any  measured  spectrum, 
the  computer  determines  the  peak-to-valloy  ratio  (the  ratio  ol  the  power 
at  the  peak  of  the  spectrum  to  that  at  the  centr  frequency)  and  the  band¬ 
width  (the  frequency  difference  measured  from  zero  Doppler  shift  to  the 
frequency  at  which  the  spectral  energy  has  fallen  by  3  dli  relative  to  the 
peak  energy).  With  these  two  values  as  abscissa  and  ordinate  of  the 
scaling  chart,  the  computer  then  determines,  from  the  skewed  grid,  3  and 

T  ;  3  is  just  the  temperature  ratio,  T  /T  ,  which  is  ini  erred  under  the 

1  2  C  3 
assumption  a  =0. 
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INCOHERENT  SCATTER  SPECTRUM  ANALYSIS  CHART  ?  =  23  cm.  PW  =  360  ps.  RG  =  *00  ^ 


•1 .  Interpretation  tor  Cr  <  1. 


Now  let  us  relax  the  restriction  we  previously  placed  on  N  and 

2  2 
thus  on  Ot  .  We  let  N  assume  its  proper  value',  and  thus  &  will  he  non- 

2 

zero,  We  require  only  Cr  ^  1 ,  This  is  a  reasonable  requirement  and  will 

2 

hold  true  throughout  most  oi  the  ionosphere.  As  cr  >1,  the  character 
o J  the  spectra  changes  drastically,  and  the  distinctive  spectral  ieatures 
d i sappear . 


The  terms  in  the  equation  describing  the  incoherent  scatter 


spectrum  are  so  arranged  that  the  spectral  shape  depends  on  o  in  the 


same  wav  that  it  depends  on  T  /-i  .  Mooreroft7  showed  ill  is  the  true 

v  e  l  r 


value  oi  the  electron-to-ion  temperature  ratio  and  p  the  fictitious  ratio 


that  would  cause  the  same  shaped  spectrum  ii  u  were  zero,  then 


T  =  (1  i  &“)& 

r 


(ii) 


This  can  be  seen  by  examining  the  equation  for  1  he  backseat  ter  spectrum. 
1  rom  l'inal  Report  2  Lp.  10,  Kq.  (lli-12'']c  we  can  write  the  denominator 
ol  the  equation  lor  the  backseutter  spectrum.  (The  numerator  does  not 


depend  on  T  or  *» ,  onlv  on  T  ,  and  thus  is  not  changed  bv  changing  T  or 

r  '  i  r 


Q.  ’  Irani  I'q.  (lll-ld)  and  assuming  singly  cliarged  ions  so  that  \’/.t)/Y. 
the  denominator,  d,  is 


d  =  [1  + 


•])"  ,  T"( 

r 


(1^) 


The  terms  indicated  by  dashes  in  parentheses  determine  the  1  requeue}  de¬ 
pendence  and  are  not  allotted  by  ii  or  Substituting  l.q.  (ll,  in 


and  t  lion  ,  collecting  terms 


Thus  we  can  see  that  the  denominator  has  the  same  lorm  lor  p 
and  a  =  0  as  it  does  lor  the  real  j  and  -1  .  This  will  enable  us,  in 

data  reduction,  to  use  only  one  spectrum  scaling  chart  (the  chart  lor 
Q  0 '  rather  than  separate  scaling  charts  depending  on  the  local  electron 
density  at  each  altitude. 

-  •; 

’1  he  magnitude  ol  the  denominator  is  increased  by  ( 1  i  Q  )  , 
thus  reducing  tin'  total  power  in  the  spectrum  by  this  same  lactor.  This 
conclusion  is  verilied  by  inspecting  the  equation  ini  the  scattering 
cross  section: 


c 

l* 

o  =  - : - - — 

[l  •  -t  T  ](1  1  q“) 
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N 


whore  is  the'  classical  radar  cross  section  oJ  a  single'  electron 

-  Li  b  ~  \  ,  Li 

(=  10  ~  nf ) .  Substituting  6, 1  -i  Q  )  lor  T  we  get 


[i  a“  (i  4  o“)g](i  t  o“)  (ir  cry (i  -i  gj 


(ig) 


which  again  shows  the  scattered  power  l educed  by  the  lactor  (  1  +  a  ) 


Now  vc  can  define  an  iterative  procedure  to  determine  ’i  ,  T  . 

e  r 

and  N,  given  'I  ,  & ,  and  N ' ,  the  raw  density  deduced  from  tin-  power  profile 
i 

under  the  assumptions  T  =  1,  g“  =  0.  The  parameters  of  interest  are 

3" 

inter lela  ted  in  tlie  following  ways: 

(l)  ri  j  and  6  are  functions  of  PV  and  I1VV  as  obtained  from 
the  spectrum  scaling  chart  for  0=0. 

(it)  Tc  -  Tr  *  Ih  =  3(1  +  02)Ti. 

(  d)  &“  =  11 .  It  2  X  10(,’i  /N. 

( -1'  N  =  [x'(l  h  &“')“(!  +  B)]/2. 


(  17 ) 
(18) 
(10) 


The  lactor  ol  2  in  hcj.  ( 1 0 )  results  fi'om  the  assumption  1  =  1,  Cr  =  0 

used  in  calculating  N ' .  The  processing  procedure  is  as  follows: 


Given ; 


Measured  parameters  N  ,  TV,  and  li\V 


Deduce : 

'1  from  t 

spectrum  scaling  chart 

■> 

i 

Initially: 

G~ 

=  0 

V 

1st  ima  t  e : 

T 

c 

from  K(|. 

(17)  \ 

\ 

A 

from  Kq. 

(19) 

Calculate : 

•> 

or 

lT'om  hcj. 

(18) 

Next ; 

Recalculate 

•> 

ric,  N,  and  ,  and  iterate 

these  three  equations  [(17),  (18),  and  (19)] 
until  a  consistent  solution  is  obtained. 

In  practice,  the  iteration  is  continued  until  both  T  and  N 

c 

change  by  less  than  1  pci'ccnt .  This  usually  takes  less  than  five 
i terat ions . 
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c. 


Assuinpt  ions ,  1  ntorpo]  at  ions  ,  ami  Kxtrapolations 


One  assumption  basic  to  the  whole  analysis  procedure  is  that  only 

"I 

one  species  oJ  ion,  namely  O  ,  is  present  in  the  ionospheric  height 
regime  over  which  we  are  probing.  This  is  probably  a  good  assumption 
over  the  height  region  where  our  radar  is  capable  of  making  spectrum 
measurements  (iiOO  to  550  km),  particularly  at  our  midlatitude  location 
at  the  present  time  near  sunspot  maximum.  The  presence  of  other  ions, 

H  I  H  H 

No  or  0  at  low  altitudes  and  He  or  at  high  altitudes,  would  not 
sign!  1  leant  ly  change  the  scattering  cross  section.  There  J'oro,  the  power 
profile  measurement  oi  raw  density  is  not  affected.  However,  significant 
amounts  ( >i(J  percent)  of  other  ions  would  markedly  affect  the  incoherent 
scatter  spectrum.  The  effect  of  multiple  ionic:  species  has  been  dis¬ 
cussed  by  MooreroM , '  lor  any  single  species,  the  shape  of  the  spectrum 
remains  essentially  unchanged ,  witli  the  width  oi  the  spectrum  being  in¬ 
versely  proportional  1o  the  square  root  oi  the  ion  mass,  lor  a  mixture 
oJ  ions,  the  spectral  shape  now  depends  on  the  relative  abundance  of  the 
species  present. 

•i 

Our  assumption  oi  only  o  ions  being  present  may  under  certain  con¬ 
ditions  lead  to  errors  in  the  deduced  temperature  at  belli  high  and  low 
altitudes.  At  the  high  altitudes,  ii  lighter  ions  were  present,  our  'I 
would  be  underestimated  and  our  i  would  be  overestimated.  At  low 
attitudes,  Ji  some  heavier  ions  were  present,  we  would  again  underestimate 

i  and  also  underestimate  'J  .  . 

l  i 

Another  assumption  used  in  the  data  processing  is  that  ho  ionospheric 
drills  or  currents  nr<-  present.  It  dri  I  Is  or  currents,  or  both,  were 
present,  they  would  cause  shills  or  asymmetries,  or  both,  in  the  measured 
spectra.  In  the  processing  procodui  •.  ,  the  spectra  are  folded  about  their 
ceb tel'  frequency,  and  the  two  sides  are  averaged  prior  to  the  determina¬ 
tion  oJ  . . .  valley  ratio  and  bandwidth.  This  is  done  to  Increase  the 
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signal- to-noisc  ratio  and  thus  improve  measurement  accuracy.  But  this 
procedure  masks  any  shifts  or  asymmetries  that  might  be  present.  At¬ 
tempts  have  been  made  to  detect  these  drifts  and  currents.  No  currents 
(spectral  asymmetries)  have  been  found.  .Small  drifts  (spectral  shifts) 
have  been  found.  Those  shifts,  when  they  are  present,  are  of  the  order 
ol  100  to  BOO  Hz  (11.5  to  B4.5  m/s).  Their  effects  on  the  temperatures 
inferred,  assuming  zero  drift,  are  very  small.  The  measured  bandwidth 
is  not  changed,  but  the  measured  peak-to-valloy  ratio  is  slightly  under¬ 
estimated,  leading  to  slight  underestimates  of  T  and  overestimates  of 
i^.  Thu  errors  in  inferred  temperature  introduced  by  drifts  of  20  m/s 
are  less  than  5  percent. 

Spectral  measurements  are  performed  only  at  eight  discrete  heights-- 
generally  200,  250,  BOO,  B50,  400,  450,  500,  and  550  lcm- -while  raw  density 
measurements  are  made  every  Jew  kilometers  between  roughly  150  and  000  km. 
In  order  to  calculate  true  densities  between  altitudes  at  which  spectra 
(temperatures)  are  available,  the  temperatures  are  linearly  interpolated 
between  the  measured  data  points.  'then,  by  using  the  interpolated  tempera¬ 
tures  and  the  raw  density,  the  true  density  can  be  computed  by  using 
).<).  (7). 

Haw  density  measurements  can  be  made  to  lower  altitudes  than  tem¬ 
perature  measurements .  Since;  received  power  is  measured  every  4,5  km 
but  the  spectral  measurement  is  made  within  a  relatively  -ong  range  gate 
(BO  to  GO  km),  it  is  common  lor  ground  clutter  to  pollute  tin;  spectral 
measurement  while  not  polluting  tin.-  raw  density  measurernei  t.  for  example, 
il  ground  el  liter  existed  out  to  a  range  oi  185  km,  the  spectrum  range 
gate  centered  at  200  l:m  would  include  some  ground  clutter,  since  it  ex¬ 
tends  i  i' on  1.82  to  218  kin .  But  the  power  prolile  measurement  would  be 
I  reo  oi  ground  clutter  lor  all  ranges  greater  than  185  km.  '1  he  lowest 
altitude  ( |(^  )  at  which  temperatures  a.  re  obtained  would  i>e  250  km  lor 
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this  ca.se,  and  some  means  of  extrapolating  the  temperatures  downward  to 
185  km  is  needed  in  order  to  obtain  true  densities  from  raw  densities 
over  the  range  185  to  250  km.  The  method  employed  makes  use  of  the 
assumption  that  the  ion  and  neutral  temperatures  arc  equal  below  R^  . 
There  is  good  theoretical  evidence'5  ” 1  °  that  this  must  be  the  case  at 
least  up  to  250  km.  The  variation  of  neutral,  and  thus  ion,  temperature 
was  assumed  to  be  of  the  form  adopted  in  the  CIKA  1965  model  atmosphere, 

T .  =  T . ( K  )  -  [ T . ( 130 ) ]  exp  [ -s( h  -  130 ) ]  ,  ( 20 ) 

i  i  L  i 

where  '1.(130  km)  is  taken  to  be  550  K,  and  the  exponent,  s,  is  given  the 
value  0.028. 

Below  R  an  assumed  form  for  T  as  a  function  of  altitude  is  also 
L  r 

needed.  The  assumption,  backed  by  evidence , 1 1  that  T  =  1  at  130  km  and 
below  is  used,  and  T  is  linearly  interpolated  between  its  value  at  R 

r  ij 

and  the  value  1  at  130  km  by  using  Kq.  (2i), 


T 


r 


1  +  ['i  (Rj  )  "  1 -°] 


~  h  -  130" 
R  -  130 

_  1 J  _ 


(21) 


Using  these  extrapolated  values 
true  density  then  found  through 


for  'i  and  T  ,  T 
r  i  e 


hq-  (7). 


can  be  computed  and  the 


At  high  altitudes  (h  >  500  km  or  so)  another  extrapolation  of  tem¬ 
peratures  is  necessary,  since  low  signal-to-noisc  ratios  limit  the  upper 
altitudes  at  which  spectral  measurements  can  be  made  before  they  limit 
the  power  profile  measurements.  The  topside  extrapolation  is  based  on  the 

assumption  that  T  is  essentially  isothermal  above  the  highest  altitude 
e 

(R  ),  at  which  spectra  are  obtained,  usually  about  500  km.  This  is  to 
be  theoretically  expected  in  the  absence  of  major  heat  sources  at  or 
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above  these  altitudes  and  because  of  tile  hi{;li  thermal  conductivity  oi 
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III  HARDWARE 

The  NASA  Project  (317  radar  used  for  the  incoherent  seat  ter  measure¬ 
ments  consists  oJ  a  number  oi'  basic  subsystems ;  the  major  subsystems  are: 
(])  transmitter  and  waveguide,  (2)  antenna  and  teed,  (3)  receivers, 

(•1J  calibration  equ  ipir.cm ,  (0)  analog- to-dig i tal  conversion,  ((3)  digital 

autooorrel a  tor,  and  ( 7 j  computer.  Several  modi i i cations  and  additions 
have  been  made  to  the  system  during  the  past  year.  This  section  describes 
the  new  equipments  and  the  new  configurations  oi  existing  equipments  im¬ 
plemented  during  the  past  year.  Table  1  reviews  the  radar  system 
parameters . 

A .  Receiving  and  Calibration  System 

In  order  to  match  the  final  IF  frequency  of  the  receiver  to  the 
digital  autocorrel ator  sampling  frequency,  several  changes  have  been  made 
in  the  receiv  "-local  oscillator  chain.  figure  11  is  a  block  diagram  oi 
the  receiver  chain.  The  signal,  alter  reception  by  the  antenna,  passes 
through  a  directional  coupler,  duplexer,  and  second  directional  coupler 
before  entering  the  parametric  amplifiers.  The  two  directional  couplers 
are  used  to  insert  calibration  signals  as  required.  (The  calibration 
signals  are  considered  further  later.)  Two  parametric  simplifiers  are 
cascaded  in  series,  resulting  in  a  110°K  system  noise  temperature,  with 
the  pa  •amps  operated  at  room  temperature  (not  cooled).  Alter  the  paramps, 
the  power  is  split  and  travels  down  two  mixer- amp]  i  J  i  er  chains,  one.'  for 
the  norm,,l  on-frequency  component  of  the  incoherent  scatter  signal  and 
one  for  the  "plasma- 1 ine"  component,  which  returns  from  the  ionosphere 
1 requency-sh i 1 1 ed  from  the  transmitted  1 requency  by  the  plasma  frequency. 
Focusing  on  the  normal  channel  (right-hand  side  oi  Figure  11),  we  see 
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I’ a  r  nine  t  e )' 

Ad  j  us  Linen  t 

Range 

No  J  in  a  1 

Opera t i on 

’J  inn.sin i  1  1  reque ney 

3280-1300  MHz 

3  200  MHz 

Transmi  t  Power  (peak) 

0-5  MW 

4  MW 

Only  Cycle-  (max.) 

r:i 

2 .  r;i 

Transmit  I’ol  a  ri  za  t  j  <)U 

RHC 

Kece  i  vl-  Pol  a  r  i  /.a  t  ion 

1,110 

Antenna  diameter 

88  i  t 

3  80  in" 

Antenna  Aperture 

Anlc-nna  Cain 

37.1  (1H 

System  Noise  Temperature 

3  3  0  "’k 

I’n  1  se  1  a  -npt  li 

3  0-500  f*s 

.'•tdO  ps 

Receiver  Video  1  iandw  i  d  l  li 

■1-200  kHz 

52  kHz 

Id  cei  vet'  1  i'  Randwi  d  t  it 

•3  0  MHz 

50  kHz 

* 


Returning  no.i  t<>  the  plasma  lnu  chain  v  1  cl  t  -  hand  side  ol  Figure  31,, 

the  power  that  was  split  oil  lol  lowing,  the  pa  ramps  enters  a  mixer.  The 

mixing  liequency  is  selected  by  the  op<  rator  so  as  to  lie  ollset  i  rum 

3320.12a  Mil/  by  the  approximate  plasma  1 requency  oi  the  1  layer.  The 

output  ol  the  mixe  r  and  tin  preamplilier  iol  lowing  the  mixe  r  is  thus  a 

band  ol  1  requeue  l  es  with  1  centered  at  30.-123  Mil/.,  where  1  is  l  he  plasma 

P  1' 

1 requency .  Alter  passing  1 rom  the  torque  tube  to  the  van,  the  signal  is 
amplil  led  and  band  1  imi  t  e  cl  to  a  1-MII/.  bandwidth  and  mixed  with  27.073  Mil/ 
to  obtain  a  signal  at  2.5  r  2  Mil/,  which  can  then  be  video-tape  recorded. 

All  local  oscillator  i requeue ies,  except  the  last  one  in  tin  plasma 
line  chain,  are  derived  1  rom  liequency  synthesizers  locked  to  a  common 
1  requency  standard.  The  clock  pulses  used  by  tin  digital  autocorrelator 
and  the  trigge  rs  lor  tile  APC’s  are  derived  1  rom  the  same  J  requency  standard. 

Provisions  lor  inserting  either  a  pul  seal  OY  calibration  signal  or  a 
pulsed  noise  source  calibration  signal  are  included  in  t  lie  receiver  sys¬ 
tem.  The  OY  calibration  is  provided  by  a  signal  generator  located  in  t lie 
van,  which  is  continuously  tunable  in  1 requency  and  continuously  adjust¬ 
able  in  amplitude.  This  signal  is  inserted  through  a  30-dIJ  directional 
coupler  immediately  l'oll owing  the  antenna  and  is  used  lor  making  Fre¬ 
quency  response1,  gain,  and  loss  measurements  throughout  the  receiver 
system.  It  is  not  used  lor  continuous  calibration  oi  t  lie  system  during 
the1  course1  oi  a  data  run,  since  the  signal  generator  3  requency  stability 
is  not  suilicie-nt  to  keep  it  within  the1  30- kHz.  band  idr  more  than  a  lew- 
tens  oi  seconds  without  constant  retuning. 

The-  pulsed  noise1  source,  inserted  through  a  20 -dll  directional  coupler 
in  Iron!  oi  the  par amps,  is  used  to  continuously  calibrate  and  monitor 
the  receiving  system  pe ri ormance  during  the  course  oi  a  data  run.  Fig¬ 
ure1  32  is  a  block  diagram  ol  the  pulsed  noise  calibration  system.  The 

,  o  \ 

system  uses  a  broadband  lie-oil  noise  source  (noise  temperature  18,000  K  J  , 


digitally  controlled  diode  switches,  and  calibrated  attenuators  to  enable 
any  of  four  selectable  amplitudes  of  wideband  pulsed  noise  to  be  injected 
into  the  receiver  chain.  With  this  system,  the  four  equivalent  noise  tem¬ 
peratures  that  can  be  used  are  70°,  115°,  17.5°,  and  8.75°K,  after  taking 
into  account  approximately  1  dll  ol  loss  introduced  by  each  diode  switch. 

Attenuation  commands  and  strobe  are  supplied  uy  the  computer.  The 
gate  supplies  pulsed  modulation  of  the  injected  noise  once  per  interpulse 
period.  The  gating  waveform  is  supplied  by  the  radar  synchronizer  and 
generated  by  counting  down  the  density  channel  ADC  trigger  pulses  to 
ensure  that  it  always  occurs  in  the  same  position  during  the  interpulse 
period.  The  countdown  circuits  are  so  arranged  that  the  gating  waveform 
always  occurs  during  the  last  sixteen  of  the  192  samples  taken  by  the  ADC . 

U .  Analog- to- Digi tal  Conversion  Equipment 

Two  ADCs  are  used  with  the  system,  one  to  digitize  the  detected 
video  signal  and  one  to  digitize  the  75-kHz  IF  signal  . 

The  density  channel  (detected  video)  data  are  digitized  by  an 
Adcom  208C  unipolar  ADC.  This  unit  is  capable  of  digitization  speeds  up 
to  100,000  samples  per  second  and  has  8-bit  resolution;  that  is,  the 
signal  is  quantized  into  256  discrete  values.  Sample  triggers  for  it 
arc  generated  in  the  Adcom  computer  interface  and  occur  at  either  30-  or 
60-gs  intervals,  selectable  by  the  operator.  The  first  sample  trigger 
occurs  at  an  operator- selected  delay  past  the  leading  edge  of  the  trans¬ 
mitted  pulse;  192  samples  arc  taken  during  each  interpulse  period. 

Figure  13  shows  the  sample  timing;  DTE  is  the  "delay  time  end,"  the  time 
at  which  the  first  sample  is  taken.  Equation  (26)  may  be  used  to  calculate 
the  apparent  radar  range  to  the  position  of  the  first  sample,  with  DN 
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FIGURE  13  DENSITY  CHANNEL  SAMPLE  TIMING 


being  tilt'  number  set  into  tlu-  thumbwheel  switches  used  to  count  clown 
the  delay: 

DTK  (km)  =  3.3  -  20  VDN  -  1,  .  (2G) 

During  the  last  sixteen  ol  the  192  sample's,  the  pulsed  noise  calibration 

signal  is  present.  The  sixteen  samples  in  lront  oi  the  calibration  pulse 

are  general  lv  used  to  estimate  the  base-line  noise  level,  P  .  Digital 

n 

numbers  produced  by  the  ADC  are  led  into  the  X1)S  930  computer  after 
rout  inn  through  the  d  i  g  i 1 a 1  autocorrelator  chassis,  which  controls  plac¬ 
ing  the  numbers  in  computer  memory. 

The  75-khz  IK  signal  is  digitized  by  a  Preston  Scientific  8500  MS 
ADC.  This  unit  is  a  9-bit  bipolar  (2’s  complement)  converter  capable  ol 
operating  at  a  maximum  speed  oi  250  khz.  The  most  significant  8  bits  of 
each  sample  are  passed  on  to  the  digital  autocorrel ator ,  which  uses  them 
in  its  autocorrelation  computation.  Triggers  for  the  Preston  ADC  are 
generated  by  the  autocorrelator.  Figure  11  shows  the  IF  channel  sample 
timing  for  a  typical  operating  setup.  Nine  bursts  of  samples  are  taken, 
eight  at  ranges  where  the  incoherent  scatter  signal  is  present  and  one 
at  a  range  where  no  signal,  only  receiver  noise,  is  present.  Eacli  burst 
consists  of  21  (or  32;  samples  spaced  by  10  us  (or  12).  The  delays 
HI,  KN ,  and  Ir  are  determined  by  one-shot  multivibrators  which  the 
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FIGURE  14  IF  CHANNEL  SAMPLE  TIMING 

operate!'  may  adjust  as  required.  The  burst  oi  samples  centered  at  |U 
is  used  by  the  autocorrelator  to  produce  tile  autocorrelation  coefficient 
lor  range  HI,  which  are  then  passed  on  to  the  XUS  9. TO  computer,  and 
similarly  ior  the  other  eight  sample  bursts, 

C .  Digi tat  Au  tocurrelator 
1  .  1  a  l  r<  ulue  t  ion 

As  part  oi  the  processing  of  incoherent  scalier  data,  it  i  s 
necessary  to  determine  the  autocorrelation  Junction  (or  equivalently, 
the  spectrum,'  oi  the  received  radar  signal  in  order  to  obtain  ionosplu  ri 
electron  and  ion  temperatures.  To  measure  tiuse  temperatures  as  a  Junc¬ 
tion  of  altitude,  the  correlation  function  must  be  measured  at  a  number 


oi  radar  ranges. 


Several  al  termillves  for  measuring  the  correlation  function  or 
spectrum  arc  possible.  A  bank  of  analog  filters  and  integrators  could 
be  implemented,  into  which  the  range- gated  received  signal  could  be  played. 
This  approach  has  several  disadvantages;  (l)  Construction  of  such  a  fil¬ 
ter  bank  is  expensive;  ( 2 )  achieving  identical  transfer  characteristics 
lor  each  filter  in  the  bank  is  dlMlcult;  (h)  achieving  identical  gains 
in  each  1 J 1  ter*  Integra  lor  channel  Is  almost  impossible;  and  ( -1 )  compensa 
ting  for  nonidentical  filler  characteristics  and  unequal  gains  in  the 
post  -  i'll  terl  ng  processing  ol  1  he  data  is  tedious  and  expensive. 

A  single  J ixed- frequency  analog  filter  could  be  used,  and  het¬ 
erodyning  technique  could  be  used  Jog  successive  analysis  ol  each  fre¬ 
quency  in  the  received  signal  through  the  same  J liter.  This  approach  h;e. 
movcj.iI  disadvantages,  slemmlng  1  rum  the  requirement  ol  range-ga t  1  ng  the 
received  signal  and  also  i  rum  the  desire  to  mlnJmi/.e  the  measurement  lime 
by  measuring  many  1  requeney  components  iiinnT  tnneously  . 

Another  possible  technique  is  to  sample  the  received  signal 
with  an  AIK'  and  then  periorm  the  auloeoi'rel  at  Jon  eilculatloii  in  n  genera) 
purpose  digital  computer  The  p'l'clpal  <11  sadvant age  ol  this  technique 
In  that  a  large  amount  ol  computer  time  is  consumed  in  the  calculation, 
making  it  Impossible  to  process  lb"  data  i  rom  more  than  a  very  lew  range 
g.ites  l  a  ici'l  time.  Thus,  ih”  data  I  loin  several  range  gates  must  Jlrsl 
be  recorded  on  digital  tap"  and  I  hen  processed  Intel.  Tills  technique  J;ns 
been  used  with  lls*  I'rojicl  bJT  radar  during  1  be  past  two  year1.,  Im j  lliree 
times  us  much  computer  time  was  reqnl.ed  to  process  the  data  as  to  record 
fll<  in . 

hurl  lip,  the  par.  f  year,  a  PJgltnl  Aui .■•eorre)  ator  was  designed 
utid  coils  l ,  iic  ted  .  Tills  device  is  essentially  a  spec  i  a) -purpose  computer 
tlint  calculates  the  mitocorro)  at  ion  Junctions  n  •  nine  ruuge  gates  in  real 
time  and  passes  the  ccrelnii-Mi  cool  J  icienls  1o  ihc  on"Jine  computer  1  or 

MU 


integration  and  recording.  The  recorded  data  can  then  be  processed 
very  quickly  (tauter  'ban  real  time),  since  the  bulk  ol  the  computations 
have  already  been  done  in  real  time  by  the  Digital  Autocorrelator . 


A  block  diagram  ol  the  incoherent  scatter  processing  system  is 
shown  in  Figure  1 A .  The  7.r>-kllz  IF  signal  is  First  bandlimited  to  Frequen¬ 
cies  between  f>0  and  loo  kHz.  and  then  sampled  at  a  100-kll/,  rate.  This 
allows  a  r>0~kllz.  hand  oi  Frequencies  to  he  analyzed,  and  "aliasing"  Folds 
the  at)-  to-10()-kllz  bund  into  the  0- 1  o-  50-klI/  region,  It  Is  Important  For 
the  ])  resampl  i  hg,  bandpass  J  11  ter  to  have  very  sharp  skirts  so  that  spurious 
1  requeney  components  it'  the  o- t  o~  50-KIIz,  band  and  at  Frequencies  above 
IDO  1.11/  do  not  appear  at  the  input  to  the  ADC,  Figure  3  0  shows  the  Fre¬ 
quency  response  oi  the  p  resampl  i  tig  J  titer  used  it;  our  system.  .Sampling, 
and  min )  og- 1  o~  di  g  1  t  a  l  conversion  oi  the  IF  signal  is  performed  by  the 
Preston  Sclent  11  le  ADC,  as  was  described  in  the  previous  section,  'liming, 
lh  the  correlator  is  generated  by  count  lug.  down  a  stable  ,'t-MII/,  clock  sig¬ 
nal  generated  by  a  I  requeney  synthes  j /.cr  ■  A  group  ol  U 1  or  'A2  samples 
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Memory  Interlace  Connector,  The  A  DCOM  ADC  samples  the  receiver's  detected 
video  signal,  which  is  used  to  generate  electron  density  profiles.  The 
sampling  oJ  the  detected  video  channel  is  asynchronous  with  the  sampling 
oi  the  1 K  channe 1 . 

The  correlator  also  receives  several  synchronizing  signals  from 
the  radar  system  and  generates  appropriate  interrupts  for  the  computer  to 
allow  the  computer's  real-time  prog, ram  to  stay  in  synchronization. 

Data  a i’e  entered  into  the  computer  through  its  Memory  Interlace 
Connector.  This  is  essentially  a  direct  pal h  to  memory  and  allows  data 
to  he  stored  in  memory  without  program  control .  A  single  computer  cycle 
Is  "stolen"  I  rum  the  running,  program  to r  each  data  word  placed  in  memory 
by  the  correlator.  The  correlator  controls  the  inputting,  oi  both  correla¬ 
tion  coeillcients  and  sampled  video  data  to  the  computer  by  genera  l  1  up. 
appropriate  computer  memory  addresses  and  control  signals, 

lor  our  appl  leal  lop,  tip-  correlator  Is  normally  so  set  up  that 
nine  bursts  (range  gates,!  ol  autocorrelation  channel  samples  are  taken  in 
each  Interpulse  period,  '1'ln-  Jlrst  burst  is  centered  at  a  nominal  range 
ol  lint)  Ian,  and  tin  next  s.'Ven  bursts  are  taken  tit  bo- km  intervals  out  to 
Dad  l.m ,  'Its  ninth  burst  Is  usually  positioned  at  a  range  greater  than 
]  ooo  l.iu,  where  no  incoherent  scalier  signal  is  present.  'J  h  I  s  allovs  the 
au toco g i a  J a l i on  Junction  ol  the  received  noise  to  be  determined  so  that 
tile  noise  can  be  siibt  meted  I  mm  the  Jirst  eight  gates  where  both  signal 
and  no  i  s< •  a  re  p n  sc  n  l  . 

As  an  example  o'  the  time  savings  made  possible  through  the  use 
o|  the  correlator,  the  el  I -line  computer  processing  v  hut  took  one 
liotii  prior  to  lla  i  mp  I  en  « ■  a  i  a  1 1  on  ol  the  correlate!'  now  lakes  about  !)  min¬ 
utes.  Tills  n  presen  l  s  an  order  ol  magnitude  savings  in  the  cost  ol  hot  It 
data  ledurtiop  personnel  and  computer  time.  Jt  u).o  greatly  reduces  the 


time  lag  between  tin  running  ol  an  experiment  and  the  availability  ol  results 


2. 


Descript ion 


The  prime  function  of  the  Digital  Autocorrelator  is  to  perform 
certain  computations  on  quantized  signals  representing  radar  returns. 

The  nature  ol  the  computations  is  such  that  they  can  be  done  by  a  general 
purpose  digital  computer,  but  only  very  inefficiently  with  respect  to  the 
time  required.  The  autocorrelator  does  these  computations  essentially  in 
real  W;e,  releasing  the  computer  for  other  calculations  at  which  it  is 
more  efficient. 

Consider  a  t ime- vary i ng  voltage,  v(t),  which  represents  a  known 

Junction  of  received  signal  amplitude  and/or  phase.  At  regular  intervals 

this  voltage  is  sampled  and  quantized  by  an  ADC.  The  quantization  is  to 

an  s-blt,  signed,  2's  compliment  binary  number,  rnd  Is  done  in  bursts  ol 

either  21  or  22  samples,  finch  lairs t  results  In  a  group  or  sequence  oJ 

numbers,  v  ,  where  l  raug.es  1  nun  1  to  !.,  and  1.  Is  either  24  or  22, 

finch  group  ol  v^'s  is  treated  as  a  separate  set  of  input  data  by  the  auto- 

coi'relator  and  it  computes  a  set  ■!  output  numbers,  A('.,  where 

X/ 

o  i  '»  I,  -  1  J 

1  N  1  •-  1.  j  .  (2V; 

1,  *»  24  or  .221 

Included  in  the  unit  are  circuitry  Jor  establishing  delays  from 
an  externally  supplied  synchron  I  z  I  ng,  signal  to  the  beg,  I  mil  ng,  ol  each  group 
oi  samples  and  lor  select  lug  I.  and  tin  sampling,  interval  (which  may  be 
e  l  the  r  1  u  o  r  1 2  .  ft  gs )  , 

'I  lie  computation  In  done  continuously,  while  samples  ol  a  group 
are  being,  accepted  and  ):  complete  within  10  ps  ol  receipt  ol  the  last 
sample.  The  transfer  ol  the  results  to  1  he  computer  takes  slightly  less 


I 
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time  than  does  the  computation,  and  a  new  group  may  commence  as  soon  as 
AC^  of  the  preceding  group  has  been  transferred. 

The  unit  is  designed  to  communicate  with  an  XDS  930  computer 
equipped  with  Multiple  Access  to  Memory  (MAM)  and  Memory  Interlace  Con¬ 
nector  (MIC;  options.  Certain  control  functions  are  accepted  from  the 
computer  as  well .  A  second  data  channel  is  built  into  the  unit,  so  that 
raw  output,  from  a  second  ADC  can  also  be  transferred  to  the  computer. 
Transfer  is  to  specified  blocks  of  memory,  and  the  two  types  of  data  are 
stored  in  separate  sets  ol  locations. 

A  1  riel  description  of  the  operation  ol  the  autocorrel a  to r 
follows,  rei eronc ) ng  the  block  diagram,  Figure  17.  Fight -bit  parallel 
data  from  the  ADC  arc  led  to  the  Input  Register,  an  8-10 t- wide- by-. 32-10 l- 
long  shill  register.  Initially,  the  register  is  cleared  to  all  zeroes. 

The  first  Fnd-ol -Conversion  (FOC;  pulse  slijlts  the  8- bit.  ADC  word  into 
the  Jirst  rank  of  the  register,  The  second  FOC  shifts  the  first  word 
down  to  tin'  second  rant.,  at  the  same  time  that  tin.'  current  word  Is  shifted 
into  the  lirst  rani.,  .subsequent  FOC  pulses  continue  to  slillt  previous 
words  down  the  register  and  place  the  new  word  into  the  first  rank.  Facli 
rani,  ol  the  Input  lleglsler  is  also  connected  to  a  corresponding  Arithmetic 
Full,  lor  which  the  register  wont  serve.'-,  as  a  mu]  t  1  pi  i  cand  , 

The  output  ol  the  AIK!  is  also  fed  to  the  Arithmetic  Control 
Fogle.  The  FOC  pulse  sets  tin  ADC  Word  Into  a  multiplier  register  In  this 
seel  ion,  otlur  circuitry  in  this  section,  triggered  by  the  FOC,  provides, 
a  sequence  ol  signals  which  cause  tin-  Arithmetic  1  nits  to  form  the  product 
ol  the  multiplier  a  no  the  respective  multiplicand  and  to  add  this  to  the 
neenmul nlrd  sum  ol  previous  products, 

Facli  ol  t  Is  Wi  identical  Arlthm'tic  Cults  contains  a  140  -  b  1 1 
accumulator  who  It  is  a  register  capable  ol  being  either  shif  ted  or  lorded 
Jii  parti)  lei.  i  In  parallel  Inputs  are  provided  by  a  30-101  parallel  adder 


r 


j 
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network,  which  forms  the  sum  oi  the  accumulator  contents  and  the  multi¬ 
plicand.  Mul tipi ication  is  done  by  an  alternating  sequence  of  adding 
into  the  accumulator  iollowed  by  shifting  the  accumulator  one  bit.  The 
adding  is  conditional  on  the  state  oi  low-order  bits  in  the  multiplier 
register.  The  result  oi  sue!  a  multiply  sequence  is  to  iorm  the  product 
ol  the  current  multiplier  and  multiplicar.fi  and  add  it  to  the  number  pre¬ 
viously  held  in  the  accumulator. 


Heiore  any  computation  starts,  all  accumulators  contain  zeroes . 

Alter  the  i i rs L  multiply  sequence,  the  J 1 rst  accumulator  contains  the 

o 

product  oi  the  it  rst  ADC  word  multiplied  by  it.se  it,  oi  .  The  remain¬ 
ing  accumulators  are  still  zero,  since  all  multiplicands  except  the  J i rst 

are  zero.  AMer  the  second  multiply  sequence,  the  1  i  rst  accumulator  con- 
2  2 

tains  v  -iv,  the  second  one  contains  v.v,,,  and  the  rest  ai".  f.ero . 
i  2  i  A 

2  2  2 

Alter  the  third  sequence,  the  accuiimi  a  tors  contain  (v  -tv  -tv), 

if  3  2  :t 


tv  v  -i  v  v  )  ,  ( v  v  )  ,  and  the  re.1 

1  2  2  1  a 


st  zeroes.  In  other  words,  arithmetic 


operations  produce  significant  results  at  first,  in  the  first,  accumulator, 
next  m  the  first  two,  then  the  i  i  rst  1  h  roe ,  etc .  Alter  the  .'{2nd ,  or  last, 
multiply  sequence,  the  accumulators  contain  the  22  AC^  which  are  now  ready 
i  o  r  t  I'a  ns  1  e  r  to  the  compu  I  e  r  . 


The  j;  in  I  Compute  signal  is  generated  at  tills  time  to  trigger  the 
Output  Control  bogie.  Here,  control  signals  are  generated  that  select  each 
accumulator  in  order  mid  do  1  he  J ol iowi ng.  first,  the  serial  output  oi  the 
accumulator  is  connected  to  the  serial  input  ol  an  output  register.  Next, 
a  series  ol  sh  i  I  t  pulses  is  applied  to  t  lie  selected  accumulator  and  to  the 
output  register,  so  as  to  tinnsler  the  accumulator  contents  to  the  output 
register,  At  the  same  lime,  the  serial  Input  to  t  lie  accumulator  is  held 
at  a  logic  zero  level,  so  that  zeroes  are  sh )  I  ted  into  the  accumulator. 

Winn  the  entire  contents  oi  'T<-  aecumu  1  ;> '  or  have  been  1  ran, si  erred  to  the 
output  register,  a  Hequcst  signal  is  sent  to  the  computer.  The  parallel 
outputs  ol  the  output  register  are  pemianently  connected,  through  line 


■Id 


drivers,  ro  the  computer  data  inputs.  When  the  computer  recognizes  the  re¬ 
quest,  it  starts  a  sequence  of  first  reading  the  address  and  then  the  data 
word.  When  the  computer  has  accepted  the  data  word,  it  sends  back  an  acknowl¬ 
edge  signal,  which  triggers  the  sequence  all  over  again  for  the  next  accumulator. 

As  has  been  pointed  out,  during  an  input- computation  cycle,  the 
accumulators  are  brought  into  significant  use  successively.  Advantage 
is  taken  of  this  to  permit  subsequent  computation  cycles  to  begin  before 
the  contents  of  all  the  accumulators  have  been  transferred  to  the  com¬ 
puter.  A  signal  is  provided  to  each  Arithmetic  Unit  (A  through  A,,) 

1  .11 

by  the  Arithmetic  Control  Logic  that  positively  inhibits  arithmetic  opera¬ 
tion  until  such  operation  can  produce  nonzero  results.  As  a  result,  some 
accumulators  can  retain  their  previously  computed  numbers,  even  while  new 
numbers  are  being  formed  in  lower- numbered  units.  The  only  limitations 
are:  (l)  that  the  subsequent  computation  cycle  not  start  until  accumulator 

J 

number  one  has  been  output  and  cleared  and  (2)  that  the  computer  accept  f 

f 

data  at  least  as  last  as  the  AIKJ  is  sampled.  The  ADC  sampling  rate  is  i 

either  100  or  S.'l  kHz  (or  one  sample  every  10  or  12  ps)  .  The  computer  i 

s 

will  accept  words  through  the  MIC  at  a  rate  ol  one  word  pci’  memory  cycle  ’ 

(•>70  kHz),  but  the  Output  Control  will  require  time  equivalent  to  three 
or  lour  cycles  to  transfer  each  word  to  the  output  register.  Tints,  an 
output  to  the  computer  will  take  place  at  intervals  oi  either  7.0  or 

H .75  ps  (d  or  5  memory  cycles).  In  either  case,  this  will  be  fast  enough  s 

f 

to  stay  ahead  ol  the  computation.  ; 

The  an  1  ocorre )  a  tor  is  packaged  in  a  5-3  /'l- i  neh- hi  gh  standard  , 

rack  mounting  unit,  A  separate  external  power  supply  provides  power  to 

the  unit.  Standard  Till  integrated  circuit  logic  is  employed,  mounted  on  ; 

1 

-10  logic  boards.  rigurc  1H  is  a  photograph  oi  the  unit  with  the  top  cover  § 

3 

% 

removed.  Figure  lit  shows  the  i  runt,  panel  of  the  autocot* relator  with  the 
operator  controls.  The  sample  rate,  number  oi  sample  hursts,  and  number 
ol  samples  pel-  burst  are  selectable  by  using  the  front  panel  toggle; 

'17 
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switches.  Potentiometers  on  the  right-hand  side  ol'  the  unit  are  used  lor 
adjusting  the  delay  to  the  first  sample,  delay  between  sample  bursts,  and 
delay  to  the  noise  burst.  Also  provided  are  a  correlator  reset  button  and 
indicator  lights  that  show  whether  the  correlator  channel  and/or  the  density 
channel  have  been  turned  on  by  the  computer.  An  operating  handbook  lor  the 
Digital  Autooo  'relator  has  been  produced,  which  describes  in  detail  the 
operation  and  maintenance  of  the  unit. 

,'i .  Programming  I  ni  onnat  ion 

.Since  the  correlator  works  in  conjunction  with  the  X1)S  9  .TO  com¬ 
puter,  certain  computer  commands  were  implemented  that  enable  interchange 
oi  information  between  correlator  and  computer. 

lour  System  ISOM  commands  are  used  to  control  operation  of  the 
correlator.  These  commands  have  the  form 

ISOM  O.'MOOX  (()  02  ,'MOOX)  , 


where  the  s i gn i J i eance  ol  X  is  as  follows: 


X  =  0 

X  r-  1 

X  =  2 

x  = 

X  =  1  through  7 


Disable  bo  til  channels. 

finable  density  channel,  disable  correlator, 

Disable  density  channel,  enable  correlator 
channel . 

finable  both  channels. 

No  el i ec t  . 


The  ISOM  is  completely  decoded  so  that  any  other  ISOM  will  have 
no  oiled.  The  computer  Start  button  disables  both  channels. 


Three  sense  output.::  are  provided  so  tlm  t  the  positions  of  the 
three  toggle  switches  can  be  determined  with  SKS  commands.  'J\vo  of  these 
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are  presently  connected.  The  command  codes  used  are: 


SKS  030001  Skip  if  32  lags,  no  skip  if  24. 

SKS  030020  Skip  if  100-kHz  sample  rate,  no  skip  if 

83.3  kHz. 

The  remaining  sense  signal  would  produce  a  skip  if  set  for  one  data 
window,  no  skip  for  eight. 

It  is  necessary  for  the  on-line  computer  program  to  stay  in 
synchronization  with  the  radar.  This  is  accomplished  through  the  use  of 
priority  interrupts.  The  interrupts  are  listed  in  Table  II  below,  by 
location  and  in  order  of  priority  (highest  first) . 


Table  II 

PRIORITY  INTERRUPT  ALLOCATION 


Interrupt 

Location 

Time  of  Occurrence 

200 

Data  Pulse  Sync 

201 

Delay  Time  End  (first  density 
channel  sample) 

202 

End  Noise  Group  (correlator) 

203 

End  each  Data  Group  (correlator) 

204 

Track  Pulse  Sync 

205 

End  ADCOM  Data  Transfer 

Correlator  output  data  are  stored  in  Memory  Locations  10000 

O 

through  10437  .  In  the  case  of  24  lags,  the  final  location  is  10327  . 

8  8 

Density  channel  data  are  stored  in  Locations  10440  through  10737  .  The 

8  8 

two  most  significant  digits  of  these  addresses  may  be  changed  by  moving 
jumper  wires  on  the  MCTL  cable  plug  module,  which  is  installed  in  the 


computer  at  Location  11B.  The  format  of  the  data  words  as  they  appear 
in  memory  is  shown  below. 


Correlator  output  word  S  -  sign  bit  (O  =  positive, 

1  =  negative) 

X  =  magnitude  bits,  2's 

SSSSSXXXXXaXXXXXXXXXXXXX  complement  form 

bit  0  4  23 

ADCOM  output  word  X  =  amplitude  bits  (unipolar) 

ooooxxxxxxxxoooooooooooo 

bit  0  4  11  23 

D.  Computer  Interfacing  Equipment 

When  the  digital  correlator  was  incorporated  into  the  radar  data 
processing  system,  a  number  of  other  changes  were  made  in  the  routing  of 
digital  data  and  conti^ol  signals  amongst  the  various  units  of  digital 
equipment.  Figure  20  is  an  interconnection  diagram  for  these  equipments, 
and  Table  III  lists  the  units  and  their  functions. 

As  part  of  this  year's  contract,  some  additional  pieces  of  computer 
equipment  were  procured  and  installed;  they  are  operating  satisfactorily. 
The  items  added  to  the  XDS  930  computer  system  include  the  following: 

(1)  Memory  Interface  Connector  option 

(2)  Multiple  Access  to  Memory  option 

(3)  Priority  Interrupt  control  chassis 

(4)  Eight  Priority  Interrupts 

(5)  Arming  feature  for  the  Priority  Interrupts. 
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Table  III 


DIGITAL  EQUIPMENTS  FUNCTION  LIST 


Un  i  t 

Functions 

ADCOM  Computer  Interface 

(1) 

Controls  XDS  930  computer  inputs  and  outputs 
through  the  PIN/POT  channel 

(2) 

Controls  pulsed  noise  source  calibrator 

(3) 

Sets  delay  to  first  sample  for  density  channel 
ADC 

(1) 

Generates  sample  triggers  for  density  channel 

AIK' 

(5) 

Controls  d igi ta 1-to-ana log  converter  and  Analok 
holding  currents 

(G) 

Provides  external  sense  line  inputs  to  computer 

(7) 

Controls  Digital  Multiplexer 

ADCOM  Digital  Multiplexer 

Gatos  21-bit  parallel  data  from  any  of  six 
selectable  data  sources  to  the  computer  inter- 
f  ace 

ADCOM  Digital-to-Analog 
Convcr ter 

Converts  digital  numbers  to  an  analog  voltage 

Analok  Circuits 

Holds  the  analog  voltages  supplied  by  tue  DAC 

Analog-to-Digi tal 

(1) 

Digitizes  detected  radar  signal 

Converters 

(2) 

Digitizes  IF  radar  signal 

Digital  Autocorrelator 

(1) 

Computes  autocorrelation  coefficients 

(2) 

Controls  transfer  of  data  to  the  XDS  930  com¬ 
puter  through  its  Memory  Interface  channel 

(3) 

Generates  interrupt  signals  for  the  computer 

d) 

Generates  status  levels  for  input  to  the  com¬ 
puter  interface 

XDS  930  Computer 

(1) 

Provides  on-line  data  processing  and  recording 
capability 

(2) 

Controls  various  portions  of  the  radar  system 

(3) 

Is  controlled  by  various  portions  of  the  radar 
system 

(!) 

Accepts  data  from  various  sources 

A  number  oi  nonstandard  computer  commands  have  been  implemented  to 
enable  the  computer  to  communicate  with  the  external  digital  equipments. 


These  commands  and  their  current  functions  are  summarized  in  Table  IV. 


Table  IV 

COMPUTER  COMMANDS 


Input  Channel 

Select  Commands 

(Follow  with  PIN  command) 

EOM  31001 

Spare 

EOM  31002 

BCD  Azimuth  (lowest  18  bits) 

BCD  units  digit  oi'  Delay  Thumbwheel  Switch  (Bits  2-5) 

EOM  31003 

BCD  elevation  (lowest  17  bits) 

BCD  tens  and  hundreds  digits  oi  Delay  Thumbwheel 

Switch  (bits  2-6) 

EOM  31004 

BCD  Time  oi  Day 

EOM  31005 

Antenna  Console  Thumbwheel  Swi tches--BCD 

EOM  3100G 

Spare 

Output  Channel 

Select  Commands 

(Follow  with  POT  command) 

EOM  32001 

Calibrate  Code  which  controls  noise  source  attenuator 

settings  (0,  1,  2,  3  give  0-,  3-,  6-,  9-dB  attenuation) 

EOM  32002 

Spare 

EOM  32003 

Spa  re 

EOM  32004 

Output  to  DAC  and  voltage  held  by  Analok  4 

EOM  32005 

Output  to  DAC  and  voltage  held  by  Analok  5 

EOM  32006 

Spare 

Table  IV  (Concluded) 


Device  Enable 

Commands 

EOM  34001 

Enable  density  channel  ADC  input  through  MIC,  disable 
autocorrelator 

EOM  34002 

Enable  autocorrelator  operation,  disable  density 

ADC  input 

EOM  34003 

Enable  both  above 

EOM  34000 

Disable  both  above 

Sense  Line 

Commands 

Op  Code 

Func  t ion 

No  Skip 

Sk  i  p 

SKS  30001 

Number  of  lags  in  each 
autocorrelator  window 

24 

32 

SKS  30002 

Sample  spacing  for 
density  channel  ADC 

30  ps 

60  ps 

SKS  30004 

Antenna  Console  Switch 
read/not  read  push¬ 
button 

DO  NOT  READ 

READ 

SKS  30010 

High-voltage  status 

OFF 

ON 

SKS  30020 

Autocorrelator  channel 
sample  spacing 

12  ps 

10  ps 

SKS  30010 


Spare 


IV  SOFTWARE 


A  large  amount  of  now  software  was  developed  during  the  Past  year 
in  order  to  make  use  of  the  new  hardware,  provide  more  operator  feedback 
during  operation,  implement  the  processing  and  analysis  procedures  de¬ 
scribed  in  Section  II,  and  speed  up  the  entire  data  recording  and  proc¬ 
essing  procedures.  These  objectives  have  been  met.  Two  major  computer 
programs  are  used:  one  that  runs  on  line  in  a  real-time  environment  and 
controls  the  data  taking  and  recording  and  one  that  is  used  off  line  to 
analyze  the  recorded  data  and  to  plot  and  print  out  the  experimentally 
determined  electron  densities  and  electron  and  ion  temperatures.  This 
sect, on  describes  these  two  programs. 

A,  On-Line  Program 
1.  Description 

This  program  is  intended  for  use  in  conjunction  with  the  data- 
tuking  facilities  of  the  XDS  930  computer  used  by  the  Project  017  radar. 

The  computer  is  presented  with  an  externally  clocked  stream  of  input  data 
(independent  of  the  computer  CPU  functioning)  consisting  of  192  amplitude 
samples  (power  profile  data)  and  24  (or  32)  lagged  correlation  coefficients 
from  each  of  nine  autocorrelation  gates  (spectral  data).  These  data  are 
input  repetitively,  once  per  radar  interpulse  period.  The  computer  is 
required  to  perform  a  running  summation  (squared  sums  for  the  amplitude 
data)  over  successive  radar  intorpulse  periods  for  each  input  number  and 
write  the  summations  in  formatted  blocks  on  magnetic  tape. 


a .  Amplitude  Samples  (Power  Profile  Data) 

The  amplitude  data  are  presented  to  the  computer  ns  n 
sequence  of  192  8-bit  unipolar  ADC  voltages  (one  sample  per  computer  word, 
Bits  4-11,  in  192  sequential  24-bit  computer  words,  from  Address  10440g 
to  10737 g) .  Each  of  these  samples  s  squared  and  added  to  on  accumulat¬ 
ing  bin  at  the  corresponding  range.  After  summation  over  some  predeter¬ 
mined  number  of  these  pulses  has  been  completed,  the  data  are  formatted 
and  output  onto  digital  tape  (if  desired),  all  the  accumulating  bins  arc 
reset  to  zero,  and  the  summation  is  repeated. 

b .  Autocorrelation  Coefficients  (Spectral  Data) 

The  autocorrelation  coefficients  are  presented  to  the 
computer  as  nine  groups  of  either  24  or  32  points  each,  read  into  sequen¬ 
tial  computer  words  (Bits  4-23)  for  a  total  of  9  X.  24  (9  X  32)  24-bit 
computer  words  per  radar  sync  pulse  (Locations  10000^  to  104378) .  Each 
of  these  groups  of  24  (32)  points  represents  a  precomputed  (by  interface 
hardware)  autocorrelation  function.  These  points  are  summed  in  9  X  24 
(9  X  32)  accumulating  bins,  and  after  a  predetermined  number  of  pulses 
t lie  data  are  output  onto  tape  (if  desired),  along  with  the  amplitude  data 
that  have  been  computed  in  a  parallel  fashion.  After  each  tape  write, 
the  autocorrelation  accumulating  bins  are  reset  to  zero,  and  the  summa¬ 
tion  is  restarted. 

c .  Control  of  the  Input  Data  Flow 

Six  priority  interrupts  are  employed  to  notify  the  program 
of  the  completion  of  data  transfers  to  the  computer  memory.  Various  sense 
lines  are  used  to  inform  the  computer  of  input  modcs--24  or  32  spectral 
samples  per  gate,  30-  or  60-gs  sample  interval  for  amplitude  data,  etc. 
(see  Table  IV) . 
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<1 .  Control  of  the  Output  Data 

The  output  from  the  program  is  in  the  form  of  digital  tope 
records,  one  record  per  summation  interval.  Each  record  contains  header 
information  ami  double  precision  summations  of  all  gates.  Tape  writing 
on  Magnetic  Tape  Unit  1  is  controlled  by  Break  Point  1  on  the  computer 
console  (set:  write  tape;  reset:  do  not  write).  The  writing  of  file 
marks  on  tape  is  controlled  by  Break  Point  4  (set:  write  file  mark; 
reset:  omit  file  mark).  The  use  of  these  switches  is  described  in  de¬ 
tail  below.  A  secondary  output  from  the  program  is  in  the  form  of  type¬ 
written  error  messages  and  visual  displays  of  the  summation  bins;  these 
are  also  described  in  detail  bo  low', 

e .  Operating  Modes 

Two  modes  of  operation  ran  be  initially  selected.  Upon 
initiation  of  the  program  at  Location  01000,  the  tvpeout --KUN(O) , 

PI/lT  (l)--is  made.  If  a  zero  is  entered  (terminated  by  a  carriage  re¬ 
turn  as  with  all  typed  input  messages),  the'  program  branches  to  the 
"run"  mode  of  operation  (described  below).  If  a  one  is  typed  in,  the 
program  enters  a  special  "check"  mode  of  operation  designed  to  chock 
individual  records  or  magnetic  tape  which  has  already  been  written,  by 
plotting  these  records  on  the*  Cal  Comp  plotter.  This  mode  is  simply  a 
convenience  for  on-the-spot  checking  of  tapes  without  loading  the  anal¬ 
ysis  program.  Either  mode  asks  the  operator  for  the  number  of  files  lie 
wishes  spaced  forward  on  MTl  with  the  message  "FILES  F1VD - ." 

l)  Plot  Mode.  After  positioning  the  tape  to  bhe  desired  file, 
the  program  types  the  message  "PEN  TO  Lll  MARGIN,  HIT  RETURN."  The  oper¬ 
ator  then  positions  the  plotter  pen  and  types  a  carriage  return.  The 

program  next  outputs  the  message'  "SAMP(O),  AC(l),  B0TH(2) - ,"  asking 

the  operator  whether  he  wishes  a  plot  of  the  192-point  amplitude  data, 


of  the  nine  autocorrelation  functions,  or  of  both.  After  accepting  a 
number  (0,  1,  or  2),  the  program  reads  one  record  from  tape,  draws  the 
desired  plots,  and  repents  the  message.  This  loop  occurs  repetitively 

until  the  program  is  restarted,  at  which  time  the  "RUN(O),  PL0T(1) - " 

message  is  typed. 


2)  Run  Mode.  Figure  21  is  a  simplified  flow  chart  of  the 
run  mode  loop.  The  run  mode  initially  executes  a  check  loop,  which  en¬ 
sures  that  MT1  is  ready  and  positioned  as  required.  Input  parameters 
are  entered  through  the  typewriter  in  response  to  typed  messages.  The 
program  then  enables  the  interface  logic  and  interrupts  and  goes  into  a 
wait  loop  seqvienco  while  processing  in  errupts.  The  wait  loop  consists 
of  a  visual  display  driver  and  a  typewriter  error  message  driver.  Errors 
are  flagged  by  the  interrupt  routines  and  typed  out  when  required.  If 
no  error  flags  are  set,  the  visual  display  is  driven  repetitively.  De¬ 
tails  of  this  operation  are  described  below. 


Console  Breakpoint  Switches 


At  the  end  of  each  sequence  of  pulses  over  which  a  summa¬ 
tion  is  performed,  the  computer  formats  the  data  for  output  onto  tape  and 
interrogates  BPl  to  see  whether  tape  output  is  desired.  If  BP1  is  set, 
the  data  are  output  onto  tap...  If  BP1  is  reset,  no  tape  write  occurs, 
unless  BPl  was  reset  during  the  current  summation  sequence.  If  that  is 
the  case,  the  buffer  is  dumped  onto  tape,  and  no  further  buffers  are  out¬ 
put  until  BPl  is  again  set.  If  BP4  is  set,  a  file  mark  is  written  fol¬ 
lowing  the  last  data  record,  and  the  block  counter  (record  numbers  within 
a  file)  is  reset  for  the  start  of  the  next  file,  when  and  Lf  it  occurs. 

An  EOF  message  is  typed,  containing  a  file  number,  time,  and  date.  If 
DPI  is  reset,  no  file  mark  is  written,  the  block  counter  is  left  set, 
and  no  message  is  typed.  This  enables  the  operator  to  maintain  block 
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continuity  within  files,  even  when  it  is  not  desired  to  have  the  blocks 
contiguous  in  time.  Unless  a  change  of  input  parameters  is  desired,  it 
should  never  by  necessary  to  restart  the  program  between  records  or  be¬ 
tween  files.  Even  when  magnetic  tape  is  not  being  written,  the  program 
remains  in  synchronization  with  the  external  interrupts;  summations  con¬ 
tinue  to  take  place;  and  program  functioning  remains  the  same  as  though 
tape  were  being  written.  When  the  first  record  in  a  new  file  has  been 
written,  a  message  is  typed  giving  the  file  number,  time,  and  date, 
thereby  generating  a  series  of  messages  defining  file  events  by  a  se¬ 
quence  number,  time,  and  date.  If  the  program  is  restarted,  the  file 
number  counters  are  not  reset  unless  MT1  is  at  the  load  point,  thereby 
maintaining  file  sequence  continuity  when  input  parameters  are  to  be 
changed.  An  EOF  should  be  written  before  restarting  the  program,  since 
the  program  initially  expects  to  begin  a  new  file  upon  startup.  If  this 
is  not  done,  a  loss  of  file  continuity  in  the  typewritten  messages  will 
result,  since  block  numbers  will  be  reset  within  a  file. 


g .  Program  Monitoring 


ent  ways: 


The  operation  of  the  program  is  monitored  in  two  differ- 


(1)  Foreground  and  background  typewriter  messages  that 
allow  the  operator  to  determine  the  time  and  sequence 
of  events  recorded  on  tape,  as  well  as  recording  the 
occurrence  of  possible  timing  errors .  Foreground 
messages  also  allow  for  the  input  of  program  param¬ 
eters  during  the  setup  phase. 

(2)  On-line  visual  monitoring  of  the  data-taking/ 
computation  process  by  means  of  ten  selectable  oscil¬ 
loscope  displays  (the  192-point  power  profile  data 
and  the  nine  autocorrelation  gates). 
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h.  Typewriter  Messages 


Several  types  of  messages  are  typed  out  during  the  running 
of  the  program.  Only  the  first  type  requires  an  operator  response.  The 
various  message  types  are  listed  below. 

(l)  Messages  controlling  the  initial  setup  of  the  run 

mode.  Certain  of  these  messages  inform  the  operator 
of  error  conditions  in  the  initial  setup  for  running: 

MT1  556  (set  MT1  to  556  bpi ) 

MT1  FPT  (no  file  ring  on  MTl) 

MT1  NOT  RDY  (MTl  not  ready) 

WRONG  AC  SAMPLE  SIZE  (program  size  does  not 
match  sample  size  of  autocorrelation  gates 
(24  or  32)  as  controlled  by  autocorrelator 
swi tch  setting ) . 

All  these  messages  may  be  cleared  by  hitting  a  car¬ 
riage  return  when  the  condition  has  been  corrected. 
The  program  repeats  the  typeouts  until  it  senses  a 
completely  correct  condition. 

Another  set  of  messages  controlling  the  initial 
setup  requests  the  operator  to  i  lput  system  param¬ 
eters,  some  of  which  are  used  tc  control  program 
operation  and  others  of  which  are  simply  output  onto 
tape  for  use  by  the  analysis  program.  These  messages 
are: 


FILES  FWD - 

SYSTEM  PRF ( 150) - 
POWER  (KW)  — 

PULSE  WIDTH( MICROSECS )  — 
CAL  IB  CODE  — 

BANDWIDTH  (KCS)  — 

AC  START  RANGE (KMS) - 

AC  NOISE  RANGE (KMS)- 


X 

i 

M 


i 


i 

:* 

1 
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date(mo/dy/yh)- * 

NO  PULSES(PRF=75) — 


Responses  to  these  messages  (with  the  exception  of 
the  date)  should  consist  of  numeric  input  only  (0-9) 
or  a  minus  sign;  they  require  termination  by  a  car¬ 
riage  return.  Any  other  character  input  results  in 
a  repetition  of  the  message.  The  date  type-in  is  a 
straight  BCD  inpu  c/  eight  characters  terminated  by 
a  carriage  return.  If  more  than  eight  characters 

are  input,  the  message  repeats.  The  FILES  FWD - 

message  accepts  •  ther  plus  or  minus  numbers.  A 
minus  one  (-2)  ped  in,  for  example,  backs  the  tape 
up  to  the  nearest  file  mark  and  then  spaces  forward 
over  the  file  r^rk,  A  plus  one  advances  the  tape 
past  the  next  file  mark.  If  the  tape  is  backed  up 
to  the  load  po.  ,  the  message 

LOAD  P0I1  .FILE  CNTRS  RESET 

is  typed,  and  1  .ntrol  is  returned  to  the  FILES  FWD - 

message.  T  ie  :.'ST  .ml  PRF(l50)-  message  accepts  only- 
even  number..!  ;  i'ers  to  the  total  PRF ,  including 

both  data  aid  ’a  k  pulses.  If  a  PRF  greater  than 
150  is  typed  i  ,  a  warning  message  is  typed  out: 

CAUTION:  PRF  NE  150;  IF  60-MICROSEC  SAMPLE 

GATE  IS  USED,  AT  LEAST  80  SAMPS  MUST  BE  READ 
IN  BEFORE  INT204  OCCURS—. 

See  SFLTC2  interrupt  descrip  .a  under  Program  Timing 
for  more  detail  on  tl  ^sage. 

(2)  Messages  specifying  the  beginning  and  end  of  files 
written  on  digital  tape: 

hh:mm:ss  mo/da/yr  BOF  xxx 

hh:mm:ss  mo/da/yr  EOF  xxx 

These  messages  contain  a  BCD  time  of  day  as  read  from 
the  time-code  generator,  the  BCD  date  as  input  during 
setup  (and  possibly  updated  during  running  when  the 
date  changes  at  0000  hours  Greenwich)  ,  and  a 


64 


three-digit  number  indicating  the  file  sequence 

number  within  a  time  (initially  set  by  the  FILES  FWD - 

message) . 

(2)  A  message  typed  out  during  running  whenever  BP2  is 
set,  which  gives  the  system  noise  temperature,  TN , 
as  computed  from  the  last  complete  record  over  which 
summation  occurred: 

hh:mm: ss  mo/da/yr  TN=xxx  K 

This  message  again  contains  the  date,  time  of  day, 
and  a  three-digit  temperature  in  degrees  Kelvin. 

(4)  Messages  flagging  the  occurrence  of  run-time  timing 
errors  or  of  various  fail  conditions: 


hh : mm : s  s 

mo/da/yr 

WRONG  AC  SAMPLE  SIZE 

hh : mm : s  s 

mo/da/yr 

REC  xxxx .PULSExxxxx  ,xx 

AC 18  INTS 

hh: mm: ss 

mo/da/yr 

REC  xxxx ,PULSExxxxx ,PREV 

SAMP  GATE  UNFIN 

hh: mm: ss 

mo/da/yr 

HVFAIL 

hh: mm: ss 

mo/da/yr 

REC  XXXX-MTl  NOT  RDY 

The  first  message  can  occur  only  when  the  ac-sample- 
size  switch  on  the  interface  box  is  inadvertently 
changed  during  running.  Since  the  memory  layout  is 
different  for  the  two  possible  sample  sizes  (24  or  32)  , 
two  separate  program  tapes  are  provided.  This  condi¬ 
tion  is  therefore  treated  as  an  irrecoverable  one  and 
results  in  a  complete  program  restart.  The  file  se¬ 
quence  counters  are  maintained,  however,  and  the 
digital  tape  may  be  repositioned  to  the  nearest  file 

mark  by  typing  in  a  minus  one  to  the  FILES  FWD - 

message.  The  second  message  flags  the  occurrence  of 
more  than  eight  interrupts  on  INT203 ,  the  autocorre¬ 
lation  interrupt  for  the  first  eight  (signal)  range 
gates.  If  this  occurs,  certain  of  the  autocorrelation 
gates  may  be  processed  before  the  values  are  actually 
updated  in  the  MIC  data  input  area.  The  third  mes¬ 
sage  flags  the  occurrence  of  an  overlap  in  sample  gate 
computations.  This  is  a  recoverable  error,  provided 
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it  is  not  flagged  on  successive  pulses.  Tape  record 
numbers  within  the  digital  tape  file  are  given  for 
each  of  these  errors,  along  with  the  pulse  number  (of 
the  total  number  of  pulses  summed  in  any  given  tape 
record)  in  which  the  error  occurred.  The  fourth 
message  flags  the  occurrence  of  a  high-voltage  fail 
in  the  transmitter.  The  computer  types  out  the  fail 
and  writes  the  partially  completed  record  on  tape 
(if  tape  is  being  written).  No  file  mark  is  written, 
regardless  of  the  state  of  BP4.  For  this  error  it  is 
assumed  that  when  the  err ir  is  cleared  the  current 
file  continues  to  bo  written;  hence,  block-count  con¬ 
tinuity  is  maintained.  A  flag  is  set  in  the  partially 
summed  record ,  and  a  true  count  of  the  number  of 
pulses  summed  over  is  recorded  on  tape.  The  1A;  fail 
condition  (tape-write  inhibited)  continues  until  13P1 
is  toggled  reset-set  (or  simply  set  if  tape  write 
was  not  in  process  when  the  HV  fail  occurred).  If 
the  computer  still  senses  an  HV  fail  condition  when 
BP1  is  toggled,  the  HV  fail  message  is  repeated  and 
the  wait  loop  reentered.  All  background  functions 
and  data  processing  remain  active  during  the  HV  fail 
period.  The  fifth  message  flags  the  occurrence  of 
an  MT1  busy  condition  during  an  attempted  tape  write. 
This  condition  is  fully  recoverable,  provided  the 
tape  has  been  set  ready  before  the  next  summation 
period  ends.  The  background  error  messages  and  visual 
display  are  suppressed,  since  the  computer  hangs  up 
in  a  tape-ready  wait  loop  while  processing  input  inter¬ 
rupts,  The  message  is  not  actually  typed  out  until 
after  the  condition  has  been  corrected,  since  the 
background  loop  is  replaced  by  the  tape-ready  loop. 

The  error  is  detectable  by  the  lack  of  background 
displays. 

The  messages  are  sequenced  in  such  a  way  as  to  provide  for 
an  orderly  record  of  events  typed  in  the  sequence  in  which  they  occur 
(typewriter  speed  being  an  obstacle  to  this  if  multiple  flags  are  set). 

All  error  messages  are  stacked  to  allow  for  the  flagging  of  multiple 
errors  having  the  same  message  typeout  (up  to  five  per  message).  The  time 
of  day  typed  out  with  each  message  is  the  time  at  which  the  typewriter 
message  is  typed  out  and  not  the  exact  time  of  occurrence  of  the  error. 
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Visual  Monitoring 


i  * 

Visual  monitoring  of  the  data-taking  process  is  executed 
as  a  background  function  and  is  part  of  the  wait  loop  (along  with  the 
typewriter  error  messages)  continually  being  executed  while  the  program 
awaits  interrupts.  Ten  displays  are  possible  in  this  loop:  the  192-point 
power  profile  or  any  of  the  nine  autocorrelation  gates.  If  the  antenna- 
console  read-do-not-read  switch  is  off,  the  computer  outputs  only  the 
192-point,  a:  ilitude  data  in  the  form  of  the  current  summation.  If  the 
read-do-not-read  switch  is  on,  the  computer  interrogates  the  antenna- 
console  thumbwheel  switch.  If  the  switch  is  set  to  C001 ,C002 ,C003 . . . , 
C009,  the  computer  displays  the  corresponding  autocorrelation  gate.  At 
any  other  setting  of  the  thumbwheel  switch,  the  computer  displays  the 
192  points  of  amplitude  data. 

Figures  22  and  23  are  photographs  of  the  visual  displays. 
Figure  22  shows  ttie  averaged  received  power  as  a  function  of  range  after 
about  40  seconds  of  integration.  The  F-layer  return  is  quite  visible, 
peaking  at  about  300  km,  as  is  ground  clutter  between  110  and  160  km. 

The  70° K  noise  pulse,  which  is  injected  for  calibration  purposes,  can  be 
seen  between  950  and  1000  km.  This  display  enables  the  operator  to  con¬ 
tinuously  monitor  the  performance  of  the  radar  system  while  taking  data. 
The  other  display  (Figure  23)  that  can  be  selected  by  the  operator  is  a 
presentation  of  the  autocorrelator  output  for  each  of  the  nine  range 
gates  whore  the  autocorrelation  function  is  computed.  This  enables 
operation  of  the  digital  autocorrelator  to  be  monitored  in  real  time 
without  interruption  of  the  data  taking.  The  visual  display  is  always 
active  regardless  of  the  error  state  (lIV  fail)  of  the  data-taking,  un¬ 
less  the  computer  is  hung  up  in  a  tape-not-ready  loop.  Thus,  by  reset¬ 
ting  BP1  before  starting  a  run,  all  the  data  processing  functions  can 
be  visually  monitored,  and  data  taking  can  bo  initiated  by  simply  set¬ 
ting  BP1.  Tape  block  numbers  are  incremented  only  when  tape  is  actually 
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written,  thereby  maintaining  block-count  continuity  within  files.  The 
displays  are  all  representations  of  the  current  state  of  the  summation 
buffers  and  are  reset  along  with  these  buffers  at  the  end  of  each  summa¬ 
tion  period.  Although  tape  writing  at  low  summation  numbers  is  not  nor¬ 
mally  desirable,  the  displays  may  be  run  at  summations  of  anywhere  from 
2  on  up.  It  should  be  kept  in  mind,  however,  that  the  display  loop  is 
subject  to  interrupt  and  that  buffers  may  be  updated  in  the  middle  of 
the  display  sweep.  This  effect  is  not  noticeable  for  larger  summations 
but  is  definitely  noticeable  over  sums  of  only  a  few  pulses.  The  visual 
display  is  output  in  the  form  of  an  XY  plot  on  the  digitnl-to-analog 
converter  and  Analok  holding  circuits.  One  onalok  contains  the  signal 
to  bo  displayed  and  the  other  a  sweep  voltage  to  drive  the  X-axis  of 
the  XY  display. 

2 .  Interrupt  Processing 

This  section  describes  in  detail  the  timing  sequence  in  which 
interrupts  are  processed  in  the  run  mode  of  operation.  Names  enclosed 
between  dashes  (-SAMPGATE-)  refer  to  labels  in  the  program  listing. 

Six  interrupts  arc  processed,  occurring  as  follows: 

I NT 2 00  Data  pulse  sync  (start  of  interpulse  period) 

INT201  Start  of  192-point  sampling  gate 

INT202  End  of  Autocorrelation  Gate  9 

INT203  End(s)  of  Autocorrelation  Gates  1  through  8 

INT204  Track  pulse  sync  (middle  of  interpulse  period) 

INT205  End  of  192-point  sampling  gate. 

The  basic  time  frame  is  taken  to  be  the  13.33-ms  (PRF  of  75)  interval 
between  any  two  consecutive  data  pulse  syncs.  For  30-gs  sample  spacing, 
Interrupts  INT201  and  INT205  are  used  to  control  processing  of  the 
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192-point  sample  gate,  and  INT204  is  ignored.  For  60-gs  sample  spacing, 
Interrupts  INT2M  and  INT205  are  used,  and  INT201  is  ignored. 

Interrupt  processing  is  such  that,  at  any  given  time,  one  of 
three  basic  interrupt  states  exists,  with  a  "catch-up"  timing  transition 
between  the  second  and  third  states  (formatting  for  tape  write).  These 
states  are  as  follows: 


a. 


State  I --Ini t 1 al lzat ion 


Lei  zero  or  more  time  frames  p  .s  as  defined  by  the  count 
number  in  -SYNCSKIP-  [five  frames  skipped  at  start  of  processing  to  get 
into  sync,  then  two  frames  skipped  for  every  tape  write  (end  of  summa¬ 
tion)]: 


INT200 

BRM 

-L0ADBLK2- 

1NT2U1 

BUM 

-N0P201- 

INT202 

BRM 

-N0P202- 

INT203 

BRM 

-NOP203- 

INT20-1 

BUM 

-N0P201- 

INT203 

BRM 

-N0P205- 

-L0ADBLK2-  counts  clown  -SYNCSKIP- ,  then  loads  State  II. 

N0P201  ,N’0P202  ,N’0P203  ,NQP2f»l  ,NQP205  clear  their  respective 
interrupts  and  take  no  other  action. 

b.  State  I I --Process! ng 

Perform  the  required  summations  of  the  autocorrelation 
gates  and  squaring-ancl-summing  of  the  192-point  sample  gate: 


INT200 

BRM 

-swe- 

INT201 

BRM 

-SAMPSTRT- 

INT202 

BRM 

-ACGATE9- 

INT203 

BRM 

-ACGATE 18- 
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I NT 2 04  BRM 


-SYNC2- 


INT205  BRM  -SAMPGATE - 

-SYNC-  chocks  current  value  of  -SYNCCNT-  for  -SYNCCT-  = 
-NXPRF-  .(no.  of  sync  pulses  passed)  =  (no.  of  sync  pulses 
typed  in)[  and  bumps  -SYNCCT-  if  not  equal.  If  equal, 

-SYNC-  reads  the  current  time,  azimuth,  and  elevation  codes 
and  loads  State  III. 

-SAMPSTRT-  checks  tlie  30-  to  G0-us  interval  sense  line. 

If  GO  us,  no  action  is  taken.  If  30  as,  -SAMPSTRT-  checks 
for  completion  of  the  previous  Sample-Rate  computation 
(-SAMP! TAG-  =  0).  It  -SAMPFIAG-  =  -1,  -SAMPGATK-  (INT205) 
has  not  finished  processing  the  previous  192  samples,  and 
an  error  is  flagged. 

-ACGATK9-  checks  for  more  titan  eight  interrupts  on  INT203 
(-ACUATE18-)  and  flags  an  error  if  -ACCOUNT-  >  8;  if 
-ACCOUNT-  ^  8,  computation  of  all  remaining  autocorrela¬ 
tion  gates  is  completed.  It  also  checks  for  -SYNCCNT-  = 
-NXPRF-  and  formats  for  output. 

-ACGATK18-  computes  one  autocorrelation  gate  (24  or  32  values) 
corresponding  to  current  value  of  -ACCOUNT-  and  bumps 
-ACCOUNT-  by  1. 

-SYNC2-  checks  30-  to  GO-gs  interval  sense  line.  If  30  gs  , 
no  action  is  taken.  If  60  us,  -SYNC2-  checks  -SAMPFLAC- 
lor  completion  of  the  previous  sample  gate  and  initiates 
background  (to  the  autocorrelation  interrupts)  computation 
of  the  first  80  samples  [minimum  number  of  samples  if  sample 
gate  ends  at  next  autocorrelation  sync  pulse  = 

(192  X  GO  -  13333/2) /GO  rs  80 J. 

-SAMPGATK-  completes  (or  begins)  computation  of  the  192- 
point  sample  gate.  If  the  sample  interval  is  30  us,  all 
192  samples  must  bo  computed  before  the  next  -SAMPSTRT- 
interrupt.  If  the  sample  interval  is  60  gs ,  the  last 
192  -  80  =  112  samples  must  be  computed  before  the  next 
-SY’NC2-  interrupt.  Since  -SAMPGATE-  is  of  lower  priority 
than  -SYNC2- ,  computation  of  the  last  112  samples  cannot 
begin  until  the  first  80  have  been  completed.  Also, 

-SYNC2-  and  -SAMPSTRT-  can  interrupt  -SAMPGATE-  and  thereby 
detect  if  processing  of  the  previous  frame  has  not  been 


completed  i a  time.  Upon  completion,  -SAMPGATE-  checks 
for  -SYNCC1.T-  =  -NXPRF-  and  initiates  formatting  for  tape 
output  if  yes.  -SAMPGATE-  remains  active  into  State  III 
if  tape  formatting  is  taking  place. 


c.  State  III — End  Processing,  Start  Output 

Complete  processing  in  current  time  frame;  format  for  tape 
output  while  skipping  two  time  frames  (13.33  ms): 


INT200 

BRM 

-LOADBLKl- 

INT201 

BRM 

-N0P201- 

INT202 

BRM 

-NOP202- 

INT203 

BRM 

-N0P203- 

INT204 

BRM 

-NOP204- 

INT205 

BRM 

-SAMPGATE- 

■LOADBLKl-  ( 

same  as 

State  i)  counts  down  -SYNCSKIP-  (lets 

two  time  frames  pass)  and  loads  State  II. 

-N0P201 ,N0P202 ,N0P203 ,N0P204-  clear  their  respective  in¬ 
terrupts  without  taking  any  other  action. 

-SAMPGATE-  is  still  active  from  State  II  and  waits  for  all 
computation  to  be  complete  prior  to  formatting  for  tape 
output.  An  error  is  flagged  if  MT1  is  not  ready.  Tape 
write  is  initiated,  range  summing  bins  are  zeroed  out  for 
the  next  summation,  and  -SYNCSKIP-  is  set  to  one  so  that 
-L0ADBLK2-  skips  two  frames  before  reloading  State  II, 

This  routine  chocks  BP1  and  BP4  and  sots  any  flags  neces¬ 
sary  for  file  start-stop  messages.  At  completion, 
-SAMPGATE-  loads  BRM  -N0P205-  into  -INT205-. 


3 ,  Mi scollaneous 

The  on-line  program  was  written  in  METASYMBOL,  the  assembly 
language  for  the  XDS  930  computer.  It  uses  approximately  6700  24-bit 
words  of  core  storage,  including  both  program  and  data  storage.  Magnetic 
tapes  written  by  this  program  arc  used  as  input  to  the  off-line  analysis 


program  described  in  Section  IV-B.  The  format  of  the  output  tape  is 
shown  in  Table  V.  Each  computer  word  shown  in  the  table  requires  four 
consecutive  6-bit.  characters  on  the  output  tape.  The  length  of  the  tape 
records  depends  on  whether  the  24-lag  or  the  32-lag  autocorrelator  option 
has  been  selected,  as  indicated  by  parentheses  for  the  last  four  entries 
in  Table  V.  Recording  is  done  in  binary  (odd  parity)  at  a  packing  den¬ 
sity  of  556  bpi ,  resulting  in  record  lengths  of  6  inches  (3344  characters) 
or  7  inches  (3920  characters).  Standard  3/4-inch  gaps  separate  records 
on  the  output  tape.  For  an  on-line  integration  period  of  1  minute,  a 
2400-ft  roll  of  computer  tape  lasts  about  70  hr. 

B.  Off-Line  Analysis  Program 

We  have  seen  in  the  preceding  section  (iV-A)  that  during  an  electron 
backscattor  radar  run,  two  sets  of  data  are  preintegrated  and  recorded 
on  tape;  the  power  samples,  also  called  density  samples,  and  the  auto¬ 
correlation  coefficients.  This  section  describes  how  the  data  are  manip¬ 
ulated  in  order  to  give,  as  a  function  of  time  and  altitude,  the  electron 
density  and  the  ion  and  electron  temperatures. 

The  input  data  to  this  program,  data  which  were  integrated  over 
relatively  short  time  periods  by  the  on-line  program,  are  integrated  over 
a  longer  period  of  time.  For  each  of  these  second  integration  periods, 
the  program  computes  the  Fourier  transform  of  the  autocorrelation  coef¬ 
ficients,  thus  giving  the  spectra  at  each  altitude.  It  next  computes  a 
raw  electron  density  curve  at  the  altitudes  of  the  192  power  samples 
("raw  density"  because  it  assumes  that  the  Debye  length  is  small  compared 
to  the  operating  wavelength) ,  The  program  also  computes  the  ion  and 
electron  temperatures,  which  are  then  used  to  recalculate  the  density,  to 
greater  accuracy,  this  time  inking  into  account  the  nonzero  value  of  ot. 
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Table  V 

MAGNETIC  TAPE  FORMAT 


Word  Number 

Contents 

1 

Record  Sequence  Number 

2 

Radar  Pulse  Repetition  Frequency  (75) 

3 

On-Line  Summation  Interval  (number  of  pulses) 

4 

Time  of  Day — BCD  codes  from  Clock 

5 

Antenna  Azimuth--hundredths  of  degrees 

G 

Antenna  Elevation--hundredths  of  degrees 

7 

Peak  Transmitted  Power — kilowatts 

8 

Pulse  Width — microseconds 

9 

Calibration  Pulse  Attonuation---0  ,3  ,6  ,9  dB 

10 

Receiver  Bandwi d th — kilohertz 

11-12 

Date — BCD  codes 

13 

Delay  Numbcr~~from  interface  thumbwheel  switch 

14 

Density  Channel  Sample  Interval — 30  or  60  p,s 

15 

Range  to  First  Autoeorrclator  Gate — kilometers 

16 

Range  Interval  Between  Gatos — kilometers 

17 

Number  of  Lags  in  Each  Gate — 24  or  32 

18  ! 

Correlator  Channel  Sample  Interval  — 10  or  12  y,s 

19 

High-Voltage  Fail  Flag:  0,  No  Fail;  1,  Fail 

20 

Range  to  Ninth  Autoeorrclator  Gato--ki lometcrs 

21-236 

(21-305) 

Low-Order  Words  of  9  y  24  (9  x  32)  Autocorrelation  Coefficients 

237-452 
(309-59 G) 

High-Order  Words  of  9  X  24  (9  X  32)  Autocorrelation  Coefficients 

453-644 

(597-788) 

Low-Order  Words  of  192  Density  Channel  Summations 

645-836 

(789-980) 

High-Order  Words  of  192  Density  Channel  Summations 

Figure  24  shows  the  main  tasks  of  this  off-line  software.  The  phys¬ 
ical  assumption  and  most  of  the  mathematical  formulas  have  been  described 
in  Section  II  of  this  report.  Details  of  the  programming  procedures  fol¬ 
low,  including  descriptions  and  explanations  of  each  of  the  subroutines. 

The  computer  used  is  an  XDS  930,  with  a  core  memory  of  16K  24-bit 
words.  Approximately  13,000  words  of  storage  are  required  by  the  program. 
The  graphs  are  plotted  during  the  computer  run  on  a  CalComp  plotter. 

Most  of  the  routines  are  coded  in  FORTRAN.  Because  of  timing  con¬ 
siderations  for  the  on-line  program,  the  input  tape  is  not  in  a  FORTRAN 
format.  Therefore,  the  input  subroutines  are  written  in  Machine  Language. 

The  time  required  to  analyze  each  integration  period  is  roughly 
5-1/2  minutes.  About  4  minutes  of  this  time  are  used  to  draw  and  label 

the  final  plot,  where  the  true  densities  and  the  temperatures  are  pre¬ 

sented  . 

1.  Description  of  the  Input  Data 

The  data  are  inputted  by  three  means. 

a.  Magnetic  Input  Tape 

The  on-line  data-taking  program,  described  previously, 
generates  the  digital  tape  that  is  the  input  of  this  program  (the  tape 

format  has  been  given  in  Table  V) .  Included  on  the  tape  are: 

(1)  A  header  of  20  words 

(2)  Twenty-four  or  32  autocorrelation  coefficients  in  each 
of  nine  range  gates  (eight  gates  where  signal  is  present 
and  one  containing  noise  only;  the  eight  signal  gates 
are  usually  Ar  •-  50  km  apart,  as  shown  in  Figure  14). 

(3)  The  192  power  samples. 
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STmRT. 


Start  Integration  Period 


^INTEGRATION  PERIOCp 
FINISHED? 


SPECTRE 


POWER 


SORTIE 


End  Integration  Period 


FIGURE  24  OFF-LINE  ANALYSIS  PROGRAM  GENERAL  FLOW  CHART 
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Table  VI 

THE  VARIABLES  IN  COMMON 


Name 

Number 

of  Words 

Definition 

D 

192 

Densi ty 

samples 

NTJII 

20 

Header 

on  output  tape 

(1) 

Last  block  number  read  in 

(2) 

PRF  (Hz) 

(3)  ) 

(4)  ) 

Start  time,  BCD  characters 

(5) 

Azimuth  (degrees  ^  100) 

(6) 

Declination  (degrees  *  100) 

(7) 

Power  (kW) 

(8) 

Pulsewidth  (jis) 

(9) 

Calib  code 

(10) 

Bandwidth  (kHz) 

(ID  | 
(12)  ) 

Date  in  BCD 

(13) 

Delay  code 

(14) 

Interval  between  density  samples  (ps) 

(15) 

First  range — RQ  (km) 

(16) 

Ar-interval  between  AC  gates  (km) 

Table  VI  (Concluded) 


Name 

Number 

of  Words 

Def ini tion 

PV 

8 

( 17 )  Number  of  AC  samples  (us) 

( 18 )  Interval  between  AC  samples  (us) 

(19)  f^F2  ( MHz  *  100) 

(20)  Number  of  pulses  integrated 

Peak-to-valley  ratios 

BW 

8 

Half -power  widths  of  spectra  (kHz) 

SN 

8 

Signal-to-noi se  ratios  of  spectra 

PC 

Calibrate  level  power 

PN 

Noise  level  power 

IB 

First  index  of  density  plot 

IE 

Last  index  of  density  plot 

MORE 

6 

Additional  parameters,  written  on  output  tape 

(1)  Record  number  on  output  tape 

(2)  ) 

.  .  ,  Last  time,  BCD  characters 

(3)  ) 

(4)  Cf  *  1000,  where  Cf  =  [N(f QF2) 1/Nmax 

(5)  ) 

:  Blanks 

(6)  1 
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b .  Typewriter  Input 


For  each  integration  time,  the  number  of  pulses,  NP ,  on 

which  the  integration  is  to  be  performed  is  typed  in,  and  also  the  value 

of  the  F-layer crit  leal  frequency,  f  F2 ,  ns  measured  bv  a  nearby  ionosonde. 

o 

Optional  input:  The  number  of  sample  points  to  be  included 
in  the  calibration  pulse  and  in  the  noise  pulse  and  the  value  of  SNMAX. 
SNMAX  is  the  lowest  signal -to-noi so  ratio  that  a  spectrum  can  have  if  it 
is  to  be  used  to  infer  temperatures.  These  numbers  are  asked  for  only 
if  Sense  Switch  1  is  set. 

c .  Spectrum  Scaling  Chart 

This  chart  (Figures  7  through  10)  is  on  paper  tape  and  is 
read  in  .iust  once  during  a  computer  run.  The  ion  temperature  and  the 
temperature  ratio  are  derived  from  this  chart,  as  has  been  seen  in 
Section  11-15. 

For  more  detail  on  these  various  inputs,  refer  to  the 
descriptions  of  each  subroutine  which  follow.  The  routines  are  described 
in  their  calling  order  (see  the  general  flow  chart,  Figure  21) •  Most  of 
the  variables  are  in  common,  and  are  list ed  with  their  definitions  in 
Table  VI. 
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2. 


MAIN  Program 


Subroutines  colled:  ENTRE 

POWER 

MAD 

Descrip  t  i  on .  To  moke  the  program  more  versatile,  the  MOVIN' 
program  calls  the  major  routines,  instead  of  having  them  linked  one  to 
the  other.  Its  flow  is  o  simple  loop: 


Add  1  to  the  output  record  count 

Call  ENTRE 

Call  POWER 

Coll  MAD 

J 


When  the  computer  run  has  ended,  the  operator  branches  to  an 
instruction  which  writes  an  end  of  file  on  the  output  tope.  Then  the 
program  comes  to  a  halt,  which,  once  cleared,  is  followed  by  o  REWIND 
of  the  output  tope. 


3 ,  Input  Routine  ENTRE 

Argumen t  s :  None 


Colled  by:  MAIN 

Routines  colled:  STARTR 
~™  CONTR 

UNPACK 

INTEGRE 

INTOF 

SATOFF 

SPECTRE 


81 


Peseri  pt  i  on .  ENTRE  is  the  ro\>tine  that  performs  the  input  of 
the  recorded  data  to  be  analyzed.  It  reads  an  appropriate  number  of 
records  from  the  input  tape,  tests  each  of  them  for  the  presence  of  a 
satellite  (in  which  case  it  drops  the  record),  and  integrates  these 
data . 

Let  us  now  follow  step  by  step  t ho  flow  chart  of  the  subroutine 
(see  Figure  25). 


a .  Tape  Ready  Test 

Thi s  is  performed  by  a  call  to  START R ,  a  machine  language 
subroutine;  the  input  tape  has  to  be  on  Un  i  t  3  and  Density  55(5. 

b .  Input  Tape  Positioning 

The  first  time  in  the  program,  one  might  want  to  start 
the  integration  on  a  record  somewhere  in  the  middle  of  the  input  tape. 

If  Sense  Switch  3  is  set,  the  time  to  position  the  tape  is  asked  for  by 
the  message:  WANTED  TIME.  After  the  wanted  time,  V.T ,  has  been  typed  in, 
the  input  tape  is  read  in  until  the  tape  time,  which  is  decoded  by  the 
Machine  Language  Routine  UNPACK,  is  equal  to  or  greater  than  WT.  The 
following  message  is  then  typed:  READ  TIME  hh  mm  sss, 

c .  First  Record  Read 

After  the  message,  NUMBER  OF  PULSES,  the  operator  types 
in  NP ,  the  number  of  pulses  on  which  he  wishes  to  perform  the  integration. 
For  a  15-minute  integration  period  and  a  PRF  of  75  per  second,  NP  =  67500, 
One  minute  is  usually  the  integration  period  used  during  the  on-line  radar 
run.  ENTRE  reads  the  first  input  recorci  into  the  buffer  IB  and  prints 
all  the  header  information  (Words  1  through  20  of  the  input  buffer).  The 
reading  routine,  called  CONTR,  is  an  entry  of  STARTR.  Then,  the  value 
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of  the  ionosonde  frequency  f  F2  is  asked  for  (see  MAD,  for  the  use  of 

o 

this  parameter). 

d .  Input  Record  Test 

A  flag  [lB(98l)T  is  set  by  CONTR ,  to  a  negative  number  if 
there  was  a  read  error,  is  set  to  a  positive  number  if  the  record  read 
was  an  END  OF  FILE  mark,  or  remains  zero  if  the  input  record  was  properly 
read . 

The  other  test  performed  here  is  the  search  for  satellites 
It  was  noticed  that  frequently  the  data  for  an  entire  15-minute  integra¬ 
tion  period  were  spoiled  by  a  satellite  passing  through  one  of  the 
radar  side  lobes  for  a  few  seconds.  This  is  why,  before  being  integrated 
the  power  profile  is  screened  by  SATOFF.  If  there  was  a  satellite,  the 
argument  IBUMP  returned  by  SATOFF  is  equal  to  the  power  sample  index  of 
the  beginning  of  the  "bun.p"  caused  by  the  reflection  on  the  satellite. 

If  IBUMP  is  not  zero,  the  data  are  not  added  to  the  accumulating  buffer, 
and  a  message  giving  the  time  and  IBUMP  is  printed. 

e.  Integrati on 

The  autocorrelation  coefficients  and  power  samples  have 
been  recorded  as  double  precision  integers,  with  the  lower-order  words 
grouped  together  first,  followed  by  the  higher-order  words  (Table  V). 
After  each  record  has  been  read,  the  autocorrelation  coefficients  and 
power  profile  sample  are  accumulated  in  Arrays  AC  and  D,  by  Machine 
Language  Subroutine  INTEGRE . 

At  the  end  of  the  integration  period  the  accumulated 
values,  AC(f,g)  and  D(j)  are  converted  into  floaing  point  variables  by 
Routine  INTOF.  However,  the  integration  is  terminated  if  an  end  of 
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file  has  been  encountered  in  the  input  tape,  or  a  high-voltage  failure 
in  the  transmitter  occurred  during  the  data  collection  period. 

In  brief,  when  ENTRE  has  finished,  Arrays  D  and  AC  contain 


(28) 


NF 

=  w  £  ACj^g)  , 


,  hero 

J  =  range  index 
i  =  pul so  indox 
i  =  lag  number 
g  =  gate  number. 

The  autocorrelation  coefficients  array  is  the  argument  of  Routine  SPECTRE 
which  is  called  at  this  point  to  compute  und  plot  the  spectra, 

\ ,  Satellite  Detection  Routine  SATQFF 

Arguments:  JH,  IRUMP,  ID  is  the  road  buffer  array; 

IJ1UMP  is  the  indox  of  the  power  sample 
whoi’c  the  satellite,  if  any,  was  detected. 

Called  by:  ENTRE 

Routines  called:  E0S1  and  E0S2 

Poser i pi  Ion,  Since  the  analysed  data  hod  often  been  spoiled 
by  an  interface  ol  a  satellite  passing  in  one  of  the  sidolobos  during 
the  radar  run,  tills  subroutine  was  written  to  screen  out  the  bud  records, 
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POWER  f JB) 


As  an  example,  Figure  26  shows  the  power,  in  decibels,  plotted  as  a  func¬ 
tion  of  the  range.  Two  satellites  were  detected  during  the  integration 
period.  The  low-altitude  satellite  spoiled  the  spectral  data,  and  the 
high-altitude  one  gave  a  wrong  noise  level.  The  importance  of  not  inte¬ 
grating  data  containing  satellite  echoes  is  apparent. 


0  100  200  300  400  500  600  700  800  900  1000 


RANGE  (km) 

FIGURE  26  POWER  PLOT  SHOWING  SATELLITE  ECHOES 

A  satellite  echo  is  a  square  pulse,  equal  in  length  to  the 
transmitted  pulse.  With  a  typical  pulse  length  of  360  ps  and  a  sample 
interval,  Ad,  of  30  ps  ,  the  number  of  samples  in  the  echo  is  T/d  =  12. 

To  screen  out  the  satellite  echoes,  we  first  test  whether  a 
detected  bump  in  the  power  profile  encompasses  eleven  or  twelve  samples. 
Since  the  satellite  echo  can  be  anywhere  in  the  power  profile,  ve  cannot 
compare  the  value  of  twelve  samples  to  the  ground  level  but  must  draw  a 
line  joining  two  power  samples  without  satellite  influence  (Line  A-B  in 
Figure  26). 


86 


Starting  at  an  altitude  of  375  km,  each  sample  D(j )  is  compared 
with  the  preceding  one;  the  difference, 


DIF  =  D(j  -  1)  -  D(j), 


(30) 


is  compared  to  SIG4 


JL  .  D(j  ~  1)  +  D(j) 
Jn 


(31) 


(N  is  the  number  of  pulses  per  input  record,  typically  4500). 

If  DIF  is  bigger  than  5a,  the  point  is  deemed  aberrant.  In 
the  case  of  a  satellite,  we  can  anticipate  two  of  these  aberrant  points, 
the  first  with  a  positive  DIF,  the  second  with  a  negative  DIF,  the  dif¬ 
ference  in  their  indexes  being  11  or  12.  This  is  what  is  done  in  the 
first  part  of  the  program  (see  Figure  27  for  the  flow  chart). 

This  test  was  found  not  to  be  sufficiently  refined,  since  re¬ 
ceived  signals  from  interfering  radars  frequently  caused  DIF  to  exceed  5a 
when  no  satellite  was  present.  Consequently,  a  further  test  was  intro¬ 
duced:  The  equation  of  the  line  AB  is  computed: 


y  =  aj  +  b 


(32) 


and  all  the  satellite  echo  samples  must  have  an  amplitude  greater  than 
aj  +  b  +  5a.  The  last  sixteen  samples  are  treated  differently  ,  because 
they  represent  the  calibration  pulse.  When  j  ^  177  (177  is  the  index  of 
the  first  calibration  sample),  the  power  D(j)  is  divided  by  the  ratio  PC/PN. 
Then  the  sample  is  submitted  to  the  above  test  procedure. 

If  the  data  do  not  satisfy  this  last  condition,  IBUMP  is  set 
to  zero  again,  and  the  remaining  samples  are  screened.  The  samples  in 
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START 


FIGURE  27  SATOFF  FLOW  CHART 


the  read-buffer  array  IB  are  not  in  a  format  suitable  to  FORTRAN.  There¬ 
fore,  a  Machine  Language  subroutine  was  written  to  convert  a  particular 
power  sample  from  the  IB  array  into  a  floating  point  number.  This  con¬ 
version  is  performed  by  Subroutines  F0S1  and  F0S2. 

5.  Spectrum  Routine  SPECTRE 

Argument:  AC(4,g) 

Called  by:  ENTRE 

Routines  called:  PTKANS 
—  pL0T 

Description .  It  has  already  been  said  (Section  III-B)  that  the 
Fourier  transform  of  the  autocorrelation  coefficients  is  the  backseat tered 
spectrum  in  which  we  are  now  interested  and  that  from  two  parameters  of 
this  spectrum — the  peak- to-valley  ratio  and  the  bandwidth--one  can  compute 
the  electron  and  ion  temperatures.  Hence,  the  three  tasks  of  SPECTRE  are: 

(1)  To  compute  the  Fourier  transforms  of  the  nine  autocorrela¬ 
tion  functions 

(2)  To  draw  the  noise  spectrum  and  the  eight  signal  spectra. 
These  plots,  in  fact,  are  here  only  lor  diagnostic  pur¬ 
poses, 

(3)  To  compute  the  PV  and  BW  of  the  eight  signal  spectra. 

The  AC  array  is  double-indexed:  SL  is  the  lag  index  (£  =  0  to  23 
or  31),  and  g  is  the  gate  index  (g  =  1  to  9) ,  with  g  =  9  corresponding  to 
the  noise  gate.  The  Fourier  cosine  transform  is  used: 

L 

S(f)  =  At  ^  AC ( f , g )  •  cos  (2rrffAt)  ,  (33) 
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where 

S  is  the  energy  at  the  frequency  f, 

At  is  the  interval  between  two  samples. 

At  limits  our  spectral  definition:  FF  =  l/(2At)  is  the  maximum 
frequency  we  can  resolve.  Numerically,  for  At  =  10  p,s,  we  have  FF  =  50  kHz. 

The  subroutine  performing  this  Fourier  transform  is  PTRANS.  It 
computes,  for  each  frequency  f, 


L-l 


S'(f)  =  AC(0fg)  +  2  £  AC(£,g)  .  cos  (  2tt£  j?At ) 

Jt=  1 


+  AC ( L , g )  •  cos  (2rtfLAt) 


(34) 


Figure  28  is  a  simplified  flow  chart  cf  this  subroutine. 


a.  Noise  Subtraction  and  Normalization 

Each  of  these  AC(4,g)  contains  both  signal  and  noise.  The 
noise  is  subtracted  and  the  coefficients  are  normalized  to  the  zero-lag 
coefficient.  The  normalization  factor,  AC(0,g)  -  AC(0,9),  is  the  total 
signal  energy  received  in  the  g  gate: 


AC(j&,g) 


AC  (  X ,  g )  -  AC(j&,9) 
AC(0,g)  -  AC(0,9) 


for  g  /  9 


(35) 


AC(£,g) 


AC(l,9) 

AC ( 0  , 9 ) 


for  g  =  9 


(36) 
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POWER 


b. 


Dandpass  Filler  Correction 


Since  the  filter  shape  is  not  a  clear-cut  battlement,  the 
spectrum  s'(f)  is  corrected  for  the  shape  of  the  noise  spectrum: 

S  (  CF ) 

noi  se 

•  *  — “ 

S 

noi  se 

This  correction  is  done  in  Subroutine  PTUANS. 

c ,  Plots 

Figure  29  is  an  example  of  the  spectra.  The  scale  is 
25  kILz/ineh  in  the  abscissa,  and  the  ordinate  is  relative  power.  The 
noise  plot  is  drawn  first,  with  f  varying  from  0  to  FF .  As  stated  pre¬ 
viously  (Section  1I-U),  the  scattered  spectra  arc  symmetric  around  zero 
Doppler  shifts,  assuming  the  effects  of  ionospheric  drifts  or  currents 
to  be  negligible.  Thus,  the  signal  is  assumed  to  be  symmetric  and  folded 
around  the  center  frequency,  CF ,  which  is  related  to  the  intermediate 
frequency,  IF,  and  to  ttic  maximum  frequency,  FF ,  by 

CF  =  2FF  -  IF  .  (38) 


S(f)  =  S'(f) 


NOISE  200  km  250  km  300  km  350  km  400  km  450  km  500  km  550  km 


FREQUENCY  —  kH; 

FIGURE  29  SAMPLE  OUTPUT  PLOTS  OF  SPECTRA 
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Typically,  CF  =  25  kHz.  The  plots  at  the*  eight  different  ranges  are  the 
"folded''  spectra;  that  is,  y  =  S(f)  +  S(2CF  -  f)  is  plotted  against  an 
abscissa  f  between  the  center  frequency  and  the  maximum  frequency. 

d.  Peak-to-Vallev  Ratio,  Bandwidth,  and  Signal- 


to-Noisc  Ratio 

During  the  computation  of  the  spectrum  S(f)  the  program 
derives  the  peak-to-valley  ratio, 

S 

max  .  , 

PV  =  -  .  /30' 

S(CF)  v  ' 

and  the  bandwidth  B\V ,  the  frequency  at  which  the  spectral  power  has 

fallen  to  S  /2.  These  two  parameters  enable  us  to  derive  the  ion.c 
max 

temperature  T  ,  and  the  temperature  ratio  5. 

To  determine  whether  our  incoherent  scatter  spectrum  is 
good'  enough  to  give  a  value  for  T  and  T  ,  the  signal-to-noi so  ratio 
(SNR)  is  used: 


......  Signal  4  Noise  „  AC(0,g) 

oMi  =  ‘  "  1  -  1  =  )  — 

Noise  AC (0,9) 


If  the  SNR  exceeds  SNMAX  (described  previously),  tlie  spectrum  is  considered 
"good"  and  is  Used  to  infer  temperatures  by  Subroutine  INTKRP, 

6.  Power  Profile  and  Raw  Densitv  Routine  POWER 


Argument : 

None 

Called  by: 

MAIN 

Subroutine  called: 

PLOT 

Output: 

Power  plot 

Descript ion .  POWER  docs  two  different  things:  It  plots  the 


received  power  (in  decibels)  as  a  function  of  the  range  (in  kilometers), 
and  it  calculates  a  raw  electron  density  using  Eq.  (3). 

Figure  2(5  is  an  example  of  the  output.  The  plot  scale  is 
2  dlJ/ineh  in  the  y  axis  and,  in  the  x  axis,  200  or  100  km/inch,  depending 
on  whether  the  interval  /d  between  the  samples  is  30  or  GO  gs. 

Figure  30  is  a  simplified  flow  chart  of  this  program.  One  sense 
switch  is  used.  If  SSI  is  set,  the  number  of  calibrate  pulse  samples, 

NC ,  and  of  noise  samples,  NN  ,  is  asked  for.  If  it  is  reset,  NC  and  NTs' 
are  left  unchanged. 

a .  Cal i brat  ion 

To  measure  the  received  power,  a  calibrated  noise  level 
is  injected  in  the  last  sixteen  density  samples.  This  calibration  pulse 
temperature  is 


T 

K 


70 


-CC/3 


(degrees  Kelvin) 


(•n) 


where  CC  is  the  calibration  code,  NUl(O).  The  exponent  CC/3  is  the  number 

of  3 -ill!  attenuators  inserted  prior  to  injection  of  the  noise  pulse  (sec 

Section  11I-A).  For  CC  =  0,  which  is  a  typical  value,  T  =  70“ ,  The 

K 

only  other  values  CC  can  have  are  3,  C> ,  or  9.  We  thus  have 


P 

K 


KT  11 
K  D 


(•12) 


where  P  is  the  power  injected  in  the  calibration  pulse,  K  is  Boltzmann's 
K 

constant,  and  15  is  the  receiver  bandwidth. 

To  determine  the  value  of  the  noise  level,  the  program 
averages  the  NN  samples  that  are  .just  before  the  N'C  calibration  samples. 


. <*•••*'-•*»  ■* . *# 


Typically,  NN  is  sixteen  samples,  but  during  the  night  periods,  when  the 
electron  density  is  very  low,  this  number  can  be  greater  for  more  pre- 
cision— 32  or  40,  During  daytime,  the  electron  density  is  sometimes  so 
high  that,  even  at  the  range  of  these  last  pulses,  some  signal  is  still 
present,  resulting  in  a  computed  noise  level  that  is  too  high. 

In  brief,  the  values  of  PC  and  PN  are  given  by 

192 

pc  =  Tc  P(J>  (,3> 

j=192--NC+l 

192-NC-l 

pk  =  ^  J2  p(j)  -  (") 

j=192-NC-NN 

b .  Power  Plot 

Along  the  Y  coordinate,  the  quantity  plotted  is 

Y(j)  =  10  log'i0  [P(j)/PN]  (dB)  .  (45) 

The  axis  is  drawn  at  Y  =  -1  dB.  The  range  along  the  X  coordinate  is 

X(j)  -  R()  +  C2(j  -  1)  (km)  ,  (40) 

where  C2  is  the  range  interval  between  samples,  and  R^  is  the  range  of 
the  first  sample.  Here  RQ  is  given  by 

R^  =  20  (D  -  1)  +  3.3  -  (km)  ,  (47) 
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where  D  is  tho  dolay  code,  t  is  tho  pulsewidth,  and  c  is  the  velocity 
of  light, 

c ,  Electron  Density 

In  this  routine  (POWER)  tho  rnw  electron  density  is  next 
computed.  I3y  combining  Eqs.  (3),  (4),  (dl),  and  (42)  wo  arrive  at  tho 
following  equation  for  raw  oloctron  donslty: 


C  ,K  70(2“C('/‘1)M  [!'(.))  - 
P  T  (PC  -  i>N) 


N  JX(.j) 


(ol/m  ) 


whore  all  tho  terms  have  previously  boon  defined.  Per  simplicity  in  the 
program,  all  the  unchanging  terms  over  an  integration  period  are  collected 
together  into  CON,  whore 


and  thus  the  raw  density  is 

D(.l)  «  CON  •  X(,l)2  •  [P(,l)  -  I»N](ol/roJ)  ,  (r>o) 

The  next  thing  to  do  is  to  compute  the  true  electron  density, 
which  is  the  purpose  of  the  program  called  MAI), 


1)7 


7. 


True  Density  Routine  MAD 


Argument:  None 

Called  by.  MAIN 

Routines  called :  I NT E RP 

SORTIE 

PLOT 

ANNOT 


Function  called:  CLI  (per forms  a  linear  interpolation) 


Desci'i  p  t.i on  .  Al  ter  going  through  Subroutine  POWF.Il ,  the 
array  I)  contains  the  raw  density.  The  values  ol  PV  and  HW  have  pre¬ 
viously  been  obtained  from  the  spectra  in  Subroutine  SPECTRE;  thus  we 
can  now  compute  a  More  Accurate  Density  (MAD). 

After  first,  initiating  the  values  of  the  parameters,  MAI)  com- 

2 

putos  at  eacli  spectrum  range  T  ,  ft,  and  v  ,  Then,  following  the  procedure 
described  in  Sections  Il-IJ-d  and  II-C,  it  computes  the  true  densities  at 

2 

each  power  sample  range  and  the  true  temperatures  corrected  for  nonzero  rj  . 
Lot  us  now  see  in  more  detail  how  these  operations  are  performed.  Fig¬ 
ure  31  is  the  flow  chart,  of  MAP , 


a,  Ini  t  i  a  1  i  z.u  t  i  on  of  Program 


Two  switches  are  used  in  this  first  purl  of  the  program: 

(1)  If  SSI  is  set,  the  value  of  SNMAX  is  asked  for; 
if  it  is  reset,  the  value  0,10  is  assumed  (if  no 
other  value  has  been  given  previously). 

(2)  If  SS2  is  set,  the  prog rum  pauses,  then  reads  from 
paper  tape  the  spectrum  scaling  chart,  This  has  to 
be  inputted  only  once,  the  first  time  the  program 
is  executed. 


OH 


DETERMINE  IB,  INDEX  OF 
THE  LOWEST  ALTITUDE 
WHERE  N  IS  COMPUTED 


DENMAX  =  0 
K  =  1 
I  -  IB 


RANGE  =  I5'  SPECTRUM  RANGE 


T,  1 

T 

eqn's 

V 

20,  21 

WRITE  TIME  &  DATE 
ON  POWER  PLOT 


T(  |  INTERP 
Tr  BETW 

Te  1  MEASUF 

OLATED 

/EEN 

?ED  VALUES 

ITERATE  TO  FIND  DO), 
THE  TRUE  DENSITY 

1 

r  (FoF2) 
W  80.fi  D 


80.6  DENMAX 


I 


FOR  I  =  IB,  IE 
D(l)  =  Cf  •  DO) 


PRINT  IB,  IE 


PRINT:  HMAX, 
foF2,  NMAX,  C, 


\  /  _ T _ 

1  1 

CALL  SORTIE 

Yes 


FIGURE  31  MAD  FLOW  CHART 


b. 


Computation  of  and  3  at  Each  "Good  Spectrum"  Altitude 


The  two  parameters  of  a  backscatter  spectrum,  the  peak- 

to-valley  ratio  and  the  bandwidth,  are  functions  of  the  temperatures  T. 

2  1 
and  T  ,  and  also  of  the  at  .  The  data  are  first  processed  by  assuming 
e 

2 

or  =0,  which  gives  a  first  value  of  T  and  of  the  temperature  ratio,  T  . 

i  r 

The  value  of  the  raw  density  is  then  used  to  correct  these  T  and  T  , 

i  r 

which,  in  turn,  are  used  to  correct  N.  After  the  raw  density  has  been 

replaced  by  this  new  N,  the  T  and  T  are  recalculated,  and  the  process 

i  r 

loops  until  consistent  values  are  found  for  N,  T  ,  and  T  .  This  itera- 

i  r 

tion  is  done  for  each  altitude  at  which  a  spectrum  has  been  computed. 

Before  being  used  to  infer  T  and  B,  a  spectrum  must 

i 

satisfy  these  four  criteria: 

(1)  9  £  BW  <  26 

(2)  1  <  PV  £  2.6 

(3)  N;  ^  0  (raw  density  at  the  altitude  of  the  con¬ 

sidered  spectrum) 

(4)  SN  >  SNMAX. 

The  index  of  the  lowest  good  spectrum  is  called  KL ,  and 
he  index  of  the  highest  one  is  KH.  For  each  K  spectrum,  KL  ^  K  ^  KH, 
Subroutine  INTERP  is  called  in,  which  uses  the  spectrum  scaling  chart 
(paper  tape  input  of  this  program)  to  interpret  PV  and  BW  in  terms  of  T 
and  0.  The  program  then  iterates  around  the  equations 


T 

r 

r  /  2. 

=  3(1+0') 

(51) 

T 

=  T.  T 

(52) 

e 

l  r 

VN 

=  ~  (1  +  a2)(l  +  a  +  Tr) 

(53) 

2 

t? 

(54) 

a 

=  14.22  x  10  T  /VN 
c* 
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until  the  temperature  T  and  the  density  converge.  At  this  point,  the 

e 

temperature  and  the  density  are  printed  out. 

When  this  loop  over  K  has  been  terminated,  KL  is  tested; 
if  it  equals  zero,  this  means  that  all  the  spectra  were  bad  and  the  den¬ 
sities  cannot  be  accurately  measured;  control  is  returned  to  the  MAIN 
program,  and  a  new  integration  is  started.  There  will  be  no  density 
plot  and  no  output  on  tape. 

c .  Computation  of  True  Density 

It  is  meaningless  to  compute  the  density  at  an  altitude 
where  ground  clutter  is  present  or  where  the  density  is  smaller  than  the 
accuracy  of  the  measurement. 

Therefore  wo  next  determine  the  altitude  extent  over  which 
to  compute  true  densities.  The  index  IB  of  the  first  . ensity  is  deter¬ 
mined  in  the  following  way.  Starting  at  the  altitude  of  the  first  good 

spectrum,  one  follows  down  the  D(j)  curve  until  D(j)  is  less  than 
10  3 

10  el/m  .  IB  is  then  set  equal  to  the  index,  j,  at  this  point. 

10 

If  D(,j )  is  never  less  than  10  ,  IB  is  sot  to  the  index  of  the  minimum 

of  the  density  curve. 

Similarly,  the  index,  IE,  of  the  highest  altitude  at 

which  the  density  is  plotted  is  determined  by  the  altitude  at  which  D(j) 

10  3  10  3 

exceeds  10  el/m  .  If  the  density  never  gets  smaller  than  10  ol/m  , 

then  IE  is  set  to  170,  equivalent  to  a  range  just  before  the  calibrate 

pulse . 

To  compute  true  densities,  T  and  T.  are  computed  at  each 

X'  1 

range  between  IB  and  IE  where  a  power  profile  sample  exists.  The  temper¬ 
atures  are  computed  differently,  depending  on  whether  the  range  is 
(l)  under  the  lowest  good  spectrum;  (2)  between  two  good  spectra;  or 
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(3)  above  the  highest  good  spectrum.  Section  II-C  describes  the  procedures 
used  for  each  of  these  three  cases. 


The  iteration  that  gives  the  true  density  from  the  raw 
density  uses  the  two  equations 


2  (i 

O'  =  11 . 22  X  10  ’  *  T  /VNP 

e 


(55) 


(1  +  o2)(l 


2 

+  O' 


+  T  ) 
r 


(56) 


where  VNP  is  the  previous  N. 

2  .  , 
To  converge  faster,  a  first  a  is  computed  from  Eq.  (55) 

and  then  an  N1  from  Eq,  (56) ,  and  the  first  VNP  is  taken  to  be  the  geo¬ 
metrical  moan  between  N '  and  Nl, 


VNP 


(57) 


When  the  density  profile  has  been  computed,  the  D  array  is 
multiplied  by  C  , 


C 

f 


N(f  F2 ) /N 
o  max 


(58) 


where 

N(f  F2)  =  (f  F2)  /80.G  (59) 

o  o 

N  =  maximum  density  computed  from  incoherent 

max 

measurement  alone, 

so  that  the  value  of  DENMAX  is  equal  to  the  one  deduced  from  1  F2, 

o 

After  the  data  have  boon  processed  by  MAD,  all  the  com¬ 
putations  for  this  integration  period  are  finished.  All  that  needs  to 

103 


!'  ilSlUlWMIlilfUiliftlMIt  flW'lllMwrt.  m  IW|>|  f  *1  WHHWIfl' !* 1 1  . . ■'-■‘tj  Ml  M Hill.  *  *  ■»  If  ■**,  r  | .  If  jfc,  .  |m:  ifia&idftyri-rH  'i  >h-  ^ 


bo  done  is  to  plot  those  values  and  to  write  them  on  the  output  tape; 
this  is  the  task  of  SORTIE,  which  is  called  at  this  point. 

8.  Output  Routine  SORTIE 

Arguments:  TI  ,'1'R 

Called  by:  MAD 

Routines  called:  PLOT 

ANNOT 

LABEL 

Description .  This  output  routine  plots  the  values  of  the 
electron  densities,  the  electron  temperatures,  and  the  ion  temperatures 
as  a  function  of  altitude.  Figure  32  is  an  example  of  this  plot. 

SORTIE  also  writes  the  computed  true  densities  and  temperatures  on  a 
magnetic  tape.  Figure  33  is  the  flow  chart  for  this  routine. 

a.  Plotting  of  the  Densities 

To  eliminati  some  noise  fluctuation,  and  since  the  trans¬ 
mitted  pulse  length  is  not  zero,  the  density  profile  is  smoothed  before 
being  plotted. 

Typically,  the  pulse  width  is  t  =  360  p,s ,  and  Ad  =  30  p,s , 
so  the  smoothing  should  be  over 

T  1 

—  X  —  =  6  samples 

2  Ad 

In  the  program  we  average  over  five  samples  and  plot  the 

averaged  density  at  the  altitude  of  the  middle  sample.  The  density  is 

11  ,  3 

plotted  on  a  linear  scale  of  2  X  10  el/m  per  inch  of  graph.  If  the 


104 


TEMPERATURE  CDEG) 


11  ,  3 

calculated  density  is  bigger  than  20  X  10  el/m  ,  the  values  are  plotted 
with  a  new  origin  shifted  10  inches  to  the  left. 


b.  Plotting  of  the  Temperatures 


The  electron  and  ion  temperatures  are  plotted  only  at  the 


altitudes  at  which  the  spectra  were  good, 


v  are  plotted  on  a  linear 


scale  of  '100  K  per  inch  of  graph.  The  measured  T  1  s  are  identified  by 

i 

an  "x"  and  the  T  1 s  by  a  small  square.  Straight  lines  connect  the  meas- 
e 

urod  temperatures.  If  one  of  the  values  is  out  of  the  scale  (T  >  3000°), 
the  pen  is  lifted,  and  the  point  is  shown  at  the  3000°  edge  of  the  graph. 


c .  Labeling  of  the  Graph 

The  graph  is  labeled  in  the  upper  center  with  the  start 
and  end  times  of  the  integration  period  and  the  date.  All  times  and  dates 
arc  GOT.  When  these  values  have  been  written,  Sense  Switch  1  is  tested. 

If  it  is  sot,  the  message 


IF  SSI  SET,  OUTPUT  TAPE  NOT  WHITTEN 


is  typed  out,  and  the  program  halts.  If,  when  the  halt  has  been  cleared, 
SSI  is  still  set,  the  labeling  is  not  finished,  and,  more  important,  the 
output  tape  is  not  written  with  the  data  corresponding  to  this  integra¬ 
tion  period. 


d .  Output  Tape 

The  tape  has  to  be  mounted  on  Unit  2,  the  density  is  556  bpi , 
and  the  tape  is  written  in  BCD  format.  This  format  has  been  adopted  so 
that  the  BCD  output  tapes  can  be  reread  by  any  other  computer  (XDS,  CDC , 
or  IBM  alike),  hopefully,' 
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l)  Positioning  of  the  Tape.  One  might  need  to  position 


the  output  tape  after  a  specific  record  in  order  to  stack  additional  data 
on  a  previously  written  output  tape.  To  do  this,  sot  both  Sense  Switches  1 
and  2.  Program  SORTIE,  after  plotting  the  temperatures  and  writing  the 
time  and  date,  tests  SSI  and  SS2.  If  both  are  sot,  it  types  the  message 

GIVE  LAST  TIME  (A4),  AND  DATE  (2A4) 


The  operator  types  in  the  last  integration  time  of  the  record  on  which  he 
wishes  to  position  the  tape.  The  tape  is  read,  record  by  record,  until; 


(1)  The  date  road  is  the  same  as  the  one  asked  for 

(2)  The  last  time  read  is  equal  to  or  bigger  than  the 
typed -in  time.  The  message  OUTPUT  TAPE  ON  hhmm  sss 
is  then  typed. 


The  next  data  will  be  written  following  this  last  read  record. 


2)  Tape  Format.  A  BCD  written  tape  has  a  physical  record 
length  of  up  to  132  characters.  One  logical  record  contains  one  or  more 
physical  records.  In  our  case,  in  one  logical  record,  the  data  correspond¬ 
ing  to  one  integration  period  are  written,  with  the  following  coding: 

WRITE  OUTPUT  TAPE  2  ,  101,  NREC,  NUII,  D,  TI ,  TR,  PC,  PN ,  IB,  IE 
l,[M0RE[j] ,J=2,6] 

101  FORMAT  [l6,  216, 2A4,  616,  2A4  ,816,  /, 

1  19  [l()E12 . 6  ,/3  ,2E12 , 6  ,  / ,  2[8E12,6,/] 

1  ,2E12 .6,16,16  ,A4 ,A3 ,  /, 316,16x1 

Table  VII  shows  the  output  tape  format. 
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Forma  l 


PRECEDING  PAOI  BLANX-NOT  FILMED 


V  SUMMARY  AND  KDCOMMKNDATIONS 

In  this  report  wo  have  explained  in  dotull  tho  uso  of  tho  Project  G17 
radar  Jor  making  Incoherent  scat  tor  measurements  of  tho  Ionosphere.  The 
duta  processing  and  analysis  methods  have  boon  fully  described,  including 
all  tho  assumptions,  interpolations,  and  extrapolations,  Those  portions 
of  the  radar  hardware  that  are  especially  important  to  tho  incoherent 
scat  ter  moasuromonts  and  those  portions  that  have  boon  developed,  modi  lied, 
or  Improved  during  the  past  year  havo  also  boon  described.  The  new 
computer  software,  in  combination  with  the  digital  autocorn  later,  lias 
resulted  In  an  ordor-ol -magn  1 1 tide  improvement  in  data  processing  lime, 
as  well  as  an  improvement  in  measurement  accuracy, 

Tho  major  improvements  incorporated  in  the  radar  system  during  the 
past  year  have  not  only  eliminated  many  inaccuracies  and  uncertainties 
that  previously  existed  hut  also  enabled  u:  1o  detect  small  ^5-  to  10- 
percent )  <  cm  I  -sysl  emit!  t  e  errors  that  were  previously  unnoticed,  The 
larger  quantities  oi  data  that  we  have  been  able  to  process  with  the 
new  system  have  also  made  us  aware  oi  the  Meed  jor  additional  computer 
programs  to  aid  in  the  i nt orpreta 1 1  on  of  the  data.  Those  two  areas 
should  be  the  object  <»j  system  improvement,  work  in  the  future.  Listed 
below  are  several  recommt uda t Ions  for  future  work  that  would  further 
improve  measurement  accuracy  and  increase  the  convenience  and  flexibility 
of  tiie  data  1  n t  erpre t a  1  i  on , 

In  comparing  the  incoherent  scatter  measurement  oi  the  1- layer 
maximum  electron  density,  NtiiiiV,  with  the  maximum  density,  N,.  as 

lull  A  1  <6 

determined  by  an  ionosonde.  we  have  Jound  small  (±0.2l»-dH)  systematic 
discrepancies  on  top  oi  a  roughly  equal  random  discrepancy.  Since  ail 
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the  incoherent  scatter  pro! lies  are  normalized  to  f  F2,  the  final 

accuracy  of  the  density  measurement  is  determined  by  the  accuracy  to 

which  1  12  can  be  scaled,  and  this  accuracy  (»s().2fi  dll)  is  the  main 
o 

cause  of  the  vundoin  discrepancies  botwoon  N  „  and  N,.  Wo  believe 

'  max 

the  systematic  d i sc ropancy  which  varies  slowly  with  time  is  due  to  a 

poor  choice  of  filter  in  the  detected  output  channel  of  the  receiver, 

A  single-polo  filter  is  used  which  causes  relatively  less  energy  to 

puss  under  conditions  ol  high  electron  temperature  (more  energy  in  the 

spectrum  wings )  than  low,  Consequently!  when  the  election  temperature 

is  high,  the  electron  density  is  underestimated.  A  redesigned  filter 

would  eliminate  tliir  problem  and  would  tend  to  make  1  he  accuracy  ol  the 

incoherent  scatter  measurement  better  Ilian  the  accuracy  to  which  i  12 

o 

can  be  sealed  j rom  the  ionosondc. 

In  thi’  area  ol  computer  soil  ware,  the  In]  lowing  items  would  be 
use! ul : 

( 1  ;  A  1 ae 1  I  It y  lor  having  dill egen t -- 1  ns  t ead  oj  the  same-- 

1  n  t  eg,  rn  1 1  on  times  lor  the  power  pro!  lie  data  and  spectral 
data,  'Jo  obtain  equivalent  accuracies,  long, eg  integration 
times  are  needed  lor  the  spectral  data  than  lor  the  power 

pro|  I  le  data, 

(2;  A  program  that  would  produce  a  motion  picture  oj  the 
denni  ly  and  t  empern  t  ur<’  pro]  lies,  I’rusontly,  these 
pro i I  leu  are  drawn  on  the  graph  plot  ter  (e,g., 

I- 1  guru  113;,  'ilils  consumes  most  ol  the  processing  time 
and  remilts  in  many  g,  mphs- - 1  ypl  ca )  ly  ah  nit  200  J  or  a 
>JH’  hour  run,  A  movie  condensing  2-1  !e  '  r  -  ol  data  into 
h  to  10  mlautus  ol  viewing,  would  be  uselu)  and  1  nl ornia 1 1  ve, 

(3;  A  I  1  ex  1  b  i  e  pa  Miia’ier  plotting,  package,  All  the  reduced 
data  are  now  available  on  one  out  put  tape,  but  no 
programs  exist  to  make  use  el  ihls  tape  in  study 
1  onospher  1 1-  parainei  i  as  .  A  desirable  program  would  b< 
one  capable  oj  reading  back  the  tape  and  plotting  any 
parameter  versus  any  other  parameter,  e.g.,  N|||;i^  versus 
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time,  N  at  altitude  h  versus  Te  at  the  same  altitude, 
scule  height  versus  height,  and  'J  r  versus  Nm,iX/N^.  yg. 

'Hits  program  would  eliminate  most  of  the  hand  plotting 
now  required  to  study  these  ionospheric  parameters. 

(  '1 )  .Segmentation  of  the  present  off-line  analysis  program 

to  allow  it  to  be  run  on  the  on-site  computer.  Currently, 
the  program  requires  almost  twice  the  core  storage 
available  in  the  on-site  computer.  31  the  program  were 
split  into  segments,  probably  it  could  be  made  to  run 
on  the  project-owned  computer.  Alternatively,  the 
computer's  memory  size  could  be  increased. 

In  summary,  the  objectives  uJ  tills  year's  program  have  been  met, 

•Si  gn  1  i  1  can  I  improvements  have  been  made  to  the  Project  017  radar  hardware 
and  software,  resulting  in  increased  use' u  1  neas  oi  the  radar  i or  studying, 
both  the  natural  and  the  disturbed  ionosphere, 
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GIOSS ARY  OF  TERMS 
KHoctivo  antenna  aperture  (m^j 

Autocorrelation  coeiTieienl  lor  lan,  •£,  and  Gate  g 
A/im’/lh  angle  oj  unl.ennu  measured  clockwise  i  rom  north  (degrees) 
HD  =  'I’l  X  Debye  lengt  h/wavelength 
1  1,22  y  3 0°T  /N  at  >.  =  0.2R  in 

V 

bandwidth  (ll/, ; 

Receiver  buudwldth  oj  correlation  channel 
Receiver  landwldlli  oj  density  cltannel 

Measured  bandwidth  ol  incoherent  scatter  spectra  1 rom  the  center 
J requeney  to  the  1 requeney  at  which  the  ) lower  Is  hall  oj  its 
peak  value  (II/.; 

Measured  el ec t ron- t o- i on  temperature  ratio  with  HD  assumed  to 
be  equal  to  zero 

Velocity  ol  light  *  It  X  10  (m/s) 

Calibration  Jailor  =  N  /N 

p  m 

12H./2  n 

.System  constant  s - - - — — — - 

ct;"0\  \ 't 1 1  y  /2  ) 
u  1/2  e  ' 

I  2  ^  l  \  1  /2 

Debye  length  =  I  K'l  e  /No  1  *>  Gtll  l  /N  )  (m) 

\  e, 1  o  /  \  e, 1  / 

Thumbwheel  switch  delay  number  to  lirst  density  channel  samjile 

Apparent  range  (on  scope)  to  lirst.  density  cltannel  sample  = 

;t,:t  -i  20. 1  IN  -  1  )Hm 


nr. 
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Ar 

At 
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JOE 
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T| 

f  K2 
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ll 

ll 

III 

Is 

I, 


I, 
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\ 


III 

C 


III 


1 


Time  interval  between  samples  for  density  channel  (either  30 
or  60  us) 

Range  interval  between  autocorrelation  range  gates  (km) 

Time  interval  between  samples  for  autocorrelation  channel 
(either  10  or  12.3  g,s ) 

-19 

Electronic  charge  =  1,6  X  10  coulomb 

Elevation  angle  of  antenna  measured  upward  item  horizon  (degrees) 

- 1 2 

Permit  tlvlty  ol  free  space  -  8,83-1  x  10  (iarads/m) 

Pt’fxluot.  oi  transmit  and  receive  system  losses 

Ordinary  wave  penetration  frequency  oi  the  12  layer  as  measured 
by  the  ionosonde 

Antenna  hall  beamw  Id  t  li  measured  i  rom  center  to  angle  in  which 
gain  is  decreased  by  3  dl!  (radians) 

Range  gate  number  ^  1  '•  g  <  9,  i or  autocorrelation  channel 
Antenna  gain  along  main  axis 
Height  above  eurlh's  surface  (m) 

Height  oj  maximum  electron  density  (m) 

-'ill 

Pol  t  xmann '  s  constant  ( Joule/deg  I. )  =  1,38  X  10 

Radar  wave  number  -  -1*VX 

Maximum  number  ol  Jags  ( 2-1  or  32 

Irtig,  number  In  an  I  ocorrej  at  J  on  computation  *•  32 

Radar  wavelength  (in; 

E  1  ee  t  ron  mass  (  Eg  J 
Jon  mass  (Eg  / 

3 

Electron  density  1  el/m  ) 

1 16 


N 


Number  of  pulses  over  which  integration  is  performed 

A  row  measured  electron  density  calculated  by  assuming 

T  =1,  oT  =  0 
r 

F2- layer  maximum  density  as  measured  by  Die  incoherent 
scatter  radar 

12-layer  maximum  density  calculated  from  i  12  =  (i  I2V/MO.6 

o  \  o  / 


Average  measured  density  channel  power  at  last  sixteen  sample 
points  (177  ,j  '•  192;,  where  injected  noise  pulse  is 
positioned  - 


±  y 

Id  t—i  j 


j  -  - 1 77 


Average  measured  density  channel  power  -  v~  ; 

Includes  both  signal  power  and  noise  power  a t  most  values  oi  j 


Injected  calibrate  noise  pulse  power  =  KT  li  (watts) 

)>  b 

Average  measured  density  channel  power  at  sixteen  sample  r.-nges 
1  mined  1  a l ely  preeediug  calibrate  pulse,  where  received  noise 
alone  (no  incoherent  scatter  signal)  Is  present  « 


Iteeeived  power  (watts) 


Feat;  transmitted  power  (watts) 


J’uise  repetition  l  requejn-y  ( II/, 


Measured  peal;- to-val  ley  latio  oj  incoherent  scatter  spectra 
it. id  a)’  range  (  m; 


I, owes l  altitude  at  wluch  measured  temperatures  are  obtained 
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HI 

HN 


SNR 

SUNK 

CT 


a 
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T 

1 


'J 

V 


T 


V 

J 


w 


Range  to  maximum  electron  density 

Uppermost  altitude  at  which  measured  temperatures  are  obtained 

Apparent  range  (on  scope)  to  center  of  first  autocorrelation 
range  gate  (km) 

Apparent  range  (on  scope)  to  center  of  autocorrelation  range 
gate  in  which  noise  alone  is  sampled 

.Signal- to- noise  ratio 

Signal-plus-noise- to-noise  ratio 


KfJoctivc*  backscatter  cross  section  of  single  electron  for 
incoherent  scatter  =  a^/^1  -t  <2^  +  T  j  j^l  +  ( m- ) 


-28 ,  2 

Rac.kscuttor  cross  section  of  single  electron  =  10  (in 


Klee tron  temperature  (degrees  K) 

Ion  temperature  (degrees  K , 

Injected  calibrate  noise  pulse  temperature  (degrees  K) 


'J  /'J  -  electron- to- ion  temperature  ratio 

e  1 

Transmitted  pulse  length  (seconds) 

tli 

Measured  density  channel  voltage  at  J  sample  point 

(1  J  1 0^  / 

Range  gate  width  (km) 
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